Soft random solids and their heterogeneous elasticity 
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Spatial heterogeneity in the elastic properties of soft random solids is examined via vulcanization 
theory. The spatial heterogeneity in the structure of soft random solids is a result of the fluctuations 
locked-in at their synthesis, which also brings heterogeneity in their elastic properties. Vulcanization 
theory studies semi-microscopic models of random-solid-forming systems, and applies replica field 
theory to deal with their quenched disorder and thermal fluctuations. The elastic deformations of 
soft random solids are argued to be described by the Goldstone sector of fluctuations contained 
in vulcanization theory, associated with a subtle form of spontaneous symmetry breaking that is 
associated with the liquid-to-random-solid transition. The resulting free energy of this Goldstone 
sector can be reinterpreted as arising from a phenomenological description of an elastic medium with 
quenched disorder. Through this comparison, we arrive at the statistics of the quenched disorder 
of the elasticity of soft random solids, in terms of residual stress and Lame-coefficient fields. In 
particular, there are large residual stresses in the equilibrium reference state, and the disorder 
correlators involving the residual stress are found to be long-ranged and governed by a universal 
parameter that also gives the mean shear modulus. 

PACS numbers: 61.43.-j,82.70.Gg,62.20.D- 
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I. INTRODUCTION 

Random solids, such as chemical gels, rubber, glasses 
and amorphous silica, are characterized by their struc- 
tural heterogeneity, which results from the randomness 
"locked-in" at the time they are synthesized. The mean 
positions of the constituent particles exhibit no long- 
range order, and every particle inhabits a unique spatial 
environment. The elasticity of random solids also inher- 
its heterogeneity from this locked-in randomness. For 
example, the Lame coefhcients and the residual stress 
vary from point to point throughout the elastic medium. 
The central goal of this paper is to develop a statistical 
characterization of random elastic media, via the mean 
values of the Lame coefficients and the residual stress 
as well as the two-point spatial correlations amongst the 
fluctuations of these quantities, which we shall name as 
the "disorder correlator" . These mean values and corre- 
lations are to be thought of as averages taken over re- 
alizations of the sample fabrication for a given set of 
fabrication parameters. We expect these characteristic 
quantities to coincide with the volume-averages of their 
single-sample counterparts. 

Our focus will be on soft random solids. These are 
network media that include chemical gels [l|, which are 
formed by the permanent random chemical bonding of 
small molecules, as well as rubber which is formed via 
the introduction of permanent random chemical cross- 
links between nearby monomers in melts or solutions of 
flexible long-chain polymers. Soft random solids are char- 
acterized by their entropic elasticity. These are media in 



which the shear modulus originates in the strong thermal 
fluctuations of the configurations of the constituent parti- 
cles and is much smaller than the bulk modulus, which is 
energetic in nature and originates in the excluded- volume 
interactions between the particles. The concept of en- 
tropic elasticity forms the basis of the classical theory 
of rubber elasticity, developed longago by Kuhn, Flory, 
Wall, Treloar and others (see Ref. g). 

As we discuss soft random solids we shall take chemi- 
cal gels as our prototype media. When the density of the 
introduced links exceeds the percolation threshold, an in- 
finite cluster of linked molecules forms, spanning the sys- 
tem, and the network acquires a thermodynamic rigidity 
with respect to shear deformations [s^] . This event is of- 
ten called the "gelation transition" or the "vulcanization 
transition" 

The geometrical or architectural aspects of the gela- 
tion/vulcanization transition can be well captured by the 
theory of percolation Q- However, to study the elastic- 
ity that emerges at the gelation/ vulcanization transition, 
and especially its heterogeneity, one needs a theory that 
incorporates not only the geometrical aspects, but also 
the equilibrium thermal fluctuations of the particle po- 
sitions and the strong, qualitative changes that they un- 
dergo at the gelation/vulcanization transition. In the 
setting of rubber elasticity, although the classical theory 
is successful in explaining the entropic nature of the shear 
rigidity of rubber, it is essentially based on a single-chain 
picture and, as such, is incapable of describing the con- 
sequences of the long scale random structure of rubbery 
media, e.g., the random spatial variations in their local 
elastic parameters and the resulting nonaffinity of their 
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local strain-response to macroscopic applied stresses. 

The general problem of heterogeneous elasticity and 
nonafSne deformations has been studied in the setting 
of flexible polymer networks and also scmi-flcxiblc 

polymer networks d, , glasses , and granular ma- 
terials [l2j . Particularly noteworthy is the recent inves- 
tigation by DiDonna and Lubensky of the general rela- 
tionship between the spatial correlations of the nonaffine 
deformations and those of the underlying quenched ran- 
dom elastic parameters [Tsj . 

The mission of the present work is to develop a statis- 
tical characterization of the heterogeneous elasticity of 
soft random solids by starting from a semi-microscopic 
model and applying a body of techniques that we shall 
call vulcanization theory to it. In particular, we aim to 
obtain the mean values and disorder correlators of elas- 
tic parameters, such as the Lame coefiicients and the 
residual stress, in terms of the parameters of the semi- 
microscopic model, such as the density of cross-links, the 
excluded-volume interactions, etc. One of our key find- 
ings is that the disorder correlator of the residual stress 
is long ranged, as arc all cross-disorder correlators be- 
tween the residual stress and the Lame coefficients. We 
also find that these disorder correlators are controlled by 
a universal scale parameter — independent of the micro- 
scopic details — that, moreover, controls the scale of the 
mean shear modulus. In addition, we characterize the 
nonaffininity of the deformations in terms of these pa- 
rameters. 

The strategy we adopt for accomplishing our goals in- 
volves a "handshaking" between two different analyti- 
cal schemes, [s^ The first scheme follows a well-trodden 
path. We begin with a semi-microscopic model, the Ran- 
domly Linked Particle Model (RLPM) fli-flT}. involv- 
ing particle coordinates and quenched random interac- 
tions between them that represent the randomness that 
is locked-in at the instant of cross-linking. In order to ac- 
count for this quenched randomness, as well as the ther- 
mal fluctuations in particle positions which are the origin 
of the entropic elasticity, we adopt the framework of vul- 
canization theory 0]. This framework includes the use 
of the replica method to eliminate the quenched random- 
ness, followed by a Hubbard-Stratonovich transformation 
to construct a field-theoretic representation in terms of 
an order parameter field — in this case, the random solidi- 
fication order parameter. We analyze this representation 
at the stationary-point level of approximation, and then 
focus on the gaplcss excitations around the stationary 
point — the Goldstone fluctuations — observing that these 
excitations can be parameterized in terms of a set of repli- 
cated shear deformation fields [13, IH, 13 • 

The second prong of our approach is less conventional. 
It begins with our introduction of a phenomenological 
model free energy of an elastic continuum, character- 
ized by a nonlocal kernel of quenched random attractions 
between mass-points. To obtain statistical information 
about this quenched random kernel, we use the replica 
method to eliminate the randomness, and obtain a pure 



model of replicas of the deformation field with couplings 
controlled by the disorder-moments of the kernel; these 
disorder-moments are then treated as unknown quanti- 
ties to be determined. This pure model has precisely the 
same structure as the Goldstone theory mentioned in the 
previous paragraph has. Thus, by comparing the two 
models we can learn the moments of the quenched ran- 
dom kernel from the (already-computed) coupling func- 
tions of the Goldstone model. Then we analyze a partic- 
ular realization of the phenomenological model having a 
fixed value of the quenched randomness. We observe that 
the natural reference state (i.e., the state of vanishing 
displacement field) of this model is not in fact an equi- 
librium state for any given realization of disorder, due 
to the random attractive interactions embodied by the 
kernel. We analyze how these attractions compete with 
the near-incompressibility of the medium to determine 
the displacement to the new, equilibrium configuration, 
which we shall term the "relaxed state'.' (This process 
can be understood in the setting of a hypothetical, in- 
stant process of preparing a sample of rubber: the cross- 
links introduce attractions and random stresses, and the 
system then undergoes relaxation, including global con- 
traction and local deformation.) We then explore shear 
deformations around this relaxed state, pass to the local 
limit, and arrive at the standard form of continuum elas- 
ticity theory, expressed in terms of the strain around the 
relaxed state, but with coefficients that are explicit func- 
tions of the quenched random kernel. Thus, using the 
information about the statistics of the quenched random 
kernel obtained via the comparison with the RLPM, we 
are able to infer statistical information about the elastic 
properties of the random elastic medium in the relaxed 
state, which is of experimental relevance. 

Why is it legitimate to identify the Goldstone theory 
arising from the microscopic model with the replica the- 
ory of the phenomenological model? The reason is that, 
within the schemes that we have chosen to analyze them, 
both models describe shear deformations not of the equi- 
librium state of the system but, rather, of the system 
immediately after cross-linking has been done but before 
any cross-linking-induced relaxation has been allowed to 
occur. The equivalence between these two schemes is 
not based solely on the equality of free energies of the 
two models; it is also based on the identity of the physi- 
cal meaning of the (replicated) deformation fields in the 
two theories, and thus the way that these fields couple to 
externally applied forces. Both schemes are descriptions 
of the elasticity of soft random solids displaying hetero- 
geneous elastic properties, one from a semi-microscopic 
viewpoint, the other invoking phenomenological param- 
eters. Thus, the equivalence of the two schemes provides 
the values of the phenomenological parameters as func- 
tions of the semi-microscopic parameters. 

Before concluding this introduction, let us emphasize 
that this work is a first attempt to derive the elastic het- 
erogeneities of vulcanized matters. To date, the prevalent 
strategy in studies of disordered systems is to assume a 
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particular structure for the quenched disorder on phe- 
nomenological grounds (such as Gaussian, short ranged, 
etc.) and explore the consequences. Assumptions about 
disorder structure are usually based on symmetry argu- 
ments and also the preference for simplicity, but oth- 
erwise lack theoretical substantiation. It is one of the 
main advantage of vulcanization theory that it can pre- 
dict some generic properties of the disordered structure 
in vulcanized matter, as is shown in the present work, 
which can be used to support and sharpen the assump- 
tions underlying more phenomenological theories. 

The classical theory of rubber elasticity, which was 
shown to be derivable from the saddle-point approxi- 
mation of vulcanization theory, is known to fail to de- 
scribe rubber elasticity in the intermediate and large de- 
formation regimes [I]. While a recent study [l^ shows 
that long wave-length thermal elastic fluctuations ac- 
count qualitatively for this failure, on general grounds, we 
expect that clastic heterogeneities should play an equally 
important role. It would therefore be interesting to ex- 
plore how the elastic heterogeneities discovered in the 
present work modify the macroscopic elasticity of rub- 
bery materials. Such a program is left for a future work. 

The outline of this paper is as follows. In Section |lll 
we analyze the semi-microscopic RLPM using the tools 
of vulcanization theory. Specifically, we use the replica 
method [20j to study a model network consisting of ran- 
domly linked particles, which exhibits a continuous phase 
transition from the liquid state to the random solid state, 
paying particular attention to the Goldstone fluctuations 
of the random solid state, which, as we have mentioned 
above, are related to the elastic shear deformations of 
the random solid state. In Section IIIIl we propose a 
nonlocal phenomenological model of a random elastic 
medium, and subsequently derive it from the RLPM, 
by identify that this phenomenological model is the low 
energy theory (i.e., it captures the Goldstone fluctua- 
tions) of the RLPM in the random solid state. Through 
this correspondence we learn information of the statis- 
tics of the quenched random nonlocal kernel. In Sec- 
tion IIVI wc study the relaxation of the phenomenolog- 
ical model to a stable state for any fixed randomness 
(i.e., any realization of disorder), due to random stresses 
and attractive interactions. We re-expand the free energy 
about this relaxed state to obtain the true elastic theory. 
This relaxed reference state is, however, still randomly 
stressed [2l|; nevertheless, the stress in this state — the 
so-called residual stress — does satisfy the condition of 
mechanical equilibrium, viz., dicrij{x) = 0. In its local 
limit, the proposed phenomenological model reproduces 
a version of Lagrangian continuum elasticity theory that 
features random Lame coefficients and residual stresses. 
In this section, we also use the phenomenological model 
to explore the related issue of elastic heterogeneity, viz., 
the nonaffine way in which the medium responds to ex- 
ternal stress. In Section |Vl we arrive at predictions for 
the statistics of the quenched random elastic parameters 
that feature in the phenomenological model in the relaxed 



state, along with the statistics of nonaffine deformations. 
Thus we provide a "first principles" account of the het- 
erogeneous elasticity of soft random solids. We conclude, 
in Section I VII with a brief summary and discussion of 
our results. 



II. SEMI-MICROSCOPIC APPROACH: THE 
RANDOMLY LINKED PARTICLE MODEL 

A. Randomly linked particle model 

The Randomly Linked Particle Model (RLPM) con- 
sists of N particles in a volume V in d dimensions. In 
order to study elasticity, including bulk deformations, V 
is allowed to fluctuate under a given pressure p. The 
positions of the particles in this fluctuating volume are 
denoted by {cj}jLi. The particles in the RLPM interact 
via two types of interactions: a repulsive interaction ve 
between all pairs of particles (cither direct or mediated 
via a solvent); and an attractive interaction between 
the pairs of particles that are chosen at random to be 
linked. We take the latter to be a soft link (as opposed to 
the usual hard constraint of vulcanization theory). Thus, 
the Hamiltonian can be written as 

JV M 
Hx= Yl ^E{C^- Cj) +J2'"Lih, - Cj J). (2.1) 

l<i<j<N e=l 

The label e, which runs from 1 to the total number of 
links M, indexes the links in a given realization of the 
quenched disorder, and specifies them via the quenched 
random variables M and {ie,je}^i- 

Wc take ve to be a strong, short-ranged repulsion, 
which serves to penalize density fluctuations and thus 
render the system nearly incompressible, as is appropri- 
ate for regular or polymeric liquids. As we shall describe 
in Section fll CI we address the interactions in Eq. (|2.ip 
by eliminating the particle coordinates in favor of collec- 
tive fields, which have the form of joint densities of the 
replicas of the particles. This continuum approach en- 
ables us to focus on the physics of random particle local- 
ization, particularly at lengthscales that are relevant for 
such localization, which (except when the density of links 
is extremely high) are long compared with the ranges of 
Ve and vl- With a focus on these longer lengthscales 
in mind, we see that it is adequate to replace the repul- 
sive interaction ve{c) by the model Dirac delta- function 
excluded- volume interaction ly^ S(c), characterized by the 
strength [2^42^ . This procedure amounts to making a 
gradient expansion in real space (or, equivalently, a wave 
vector expansion in Fourier space) of ve and retaining 
only the zeroth-order term; it gives for the strength 
the value JdcvE{c). Terms of higher order in the gra- 
dient expansion would have a non-negligible impact on 
the suppression of density fiuctuations only at length- 
scales comparable to or shorter than the range of ve, 
and fluctuation modes at such lengthscales are not the 
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ones driven via random linking to the instability associ- 
ated with random localization (and thus are not modes 
in need of stabilization via ve)- We remark that the ap- 
proximate interaction S{c) is not, in practice, singular, 
and is instead regularized via a high wave vector cut-off. 

At our coarse-grained level of description, the particles 
of the RLPM can be identified with polymers or small 
molecules, and the soft links can be identified with molec- 
ular chains that bind the molecules to one another. The 
potential for the soft links can be modeled as Gaussian 
chains 
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2a2 ■ 



(2.2) 



i.e., a harmonic attraction, or a "zero rest-length" spring, 
of lengthscale a between the two particles. In making this 
coarse-graining one is assuming that microscopic details 
(e.g., the precise locations of the cross- links on the poly- 
mers, the internal conformational degrees of freedom of 
the polymers, and the effects of entanglement) do not 
play significant roles for the long- wavelength physics. In 
part, these assumptions are justified by studying more 
detailed models, in which the conformational degrees of 
freedom of the polymers are retained 0. However, we 
should point out that the precise form oi vl is not im- 
portant for long-wavelength physics and, hence, for the 
elastic properties that wc are aiming to investigate [cf. the 
discussion at the end of Sec. Ill C] . 

From the discussion above, the RLPM is a convenient 
minimal model of soft random solids, inasmuch as it ade- 
quately captures the necessary long-wavclcngth physics. 
It can be regarded as either a model of a chemical gel, 
or as a caricature of vulcanized rubber or other soft ran- 
dom solid. The RLPM can be viewed as a simplified 
version of vulcanization theory [isl . [25| , with microscopic 
details, such as polymer chain conformations, being ig- 
nored. Nevertheless, it is able to reproduce the same 
universality class as vulcanization theory at the liquid- 
to-random-solid transition. For the study of elasticity, 
we shall consider length-scales on which the system is a 
well-defined solid (i.e., scales longer than the "localiza- 
tion length" as we shall see later in this paper). Both 
of these scales are much larger than the characteristic 
linear dimension of an individual polymer. The RLPM 
is a model very much in the spirit of lattice percolation, 
except that it naturally allows for particle motion as well 
as particle connectivity, and is therefore suitable for the 
study of continuum elasticity and other issues associated 
with the (thermal or deformational) motion of the con- 
stituent entities. 

Equation (|2.1[) is a Hamiltonian for a given realization 
of quenched disorder x = {"i-e, je}eLi; which describes 
the particular random instance of the linking of the par- 
ticles. These links are the quenched disorder of the sys- 
tem, which are specified at synthesis and do not change 
with thermal fiuctuations. This is because there is a 
wide separation between the timescalc for the linked- 
particlc system to reach thermal equilibrium and the 



much longer timescale required for the links themselves 
to break. Therefore, we treat the links as permanent. 
Later, we shall apply the replica technique [20| to aver- 
age over these permanent random links. 



B. Replica statistical mechanics of the RLPM 

For a given volume and a given realization of disorder 
X wc can write the partition function for the RLPM 
as 

. N 

^ i—l 

= Zc{V){ nA^dc, -c,J)) , (2.3) 



where Hq = ^ ^Yl!i j=\ ^('^i ~ '^j) *he excluded-volume 
interaction part of the Hamiltonian, and Zc{V) = 
Jv IliLi dci exp ( — Ho/ksT^ is the partition function of 
the liquid in the absence of any links. The issue of the 
Gibbs factorial factor, which is normally introduced to 
compensate for the overcounting of identical configura- 
tion, is a genuinely subtle one in the context of random 
solids (for a discussion, see Ref. @). However, our focus 
will be on "obscrvables" such as order parameter rather 
than on free energies, and thus the omission of the Gibbs 
factor is of no consequence. The factor 



(2.4) 



is associated with the link-induced attractive interaction 
term in the Hamiltonian. The average (• • ■)f", taken with 
respect to a Boltzmann weight involving the excluded- 
volume interaction Hamiltonian Hq, is defined as 



Ho 



ZciV) 



N 



(2.5) 



The corresponding Hclmholtz free energy is then given 

by 



F^iV)^-kBTlnZ^{V). 



(2.6) 



To perform the average of the free energy over the 
quenched disorder, we shall need to choose a probabil- 
ity distribution that assigns a sensible statistical weight 
Vdie, je}tLi) to each possible realization of the total 
number M and location {ie,je}tLi of the links. Follow- 
ing an elegant strategy due to Deam and Edwards [23j . 
we assume a version of the normalized link distribution 
as follows: 



nx) 



(i^)"^x(^o) 



MlZi 



(2.7) 
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where rf^ is a parameter that controls the mean total 
number of links. We assume that the preparation state 
(i.e., the state in which the links are going to be intro- 
duced) is in a given volume Vq- The Z^{Vo) factor is actu- 
ally the partition function, as given in Eq. (|2.3p , and can 
be regarded as probing the equilibrium correlations of the 
underlying unlinked liquid. The factor Aq*^-* = (27ra^)''^^ 
is actually the p — value of the Fourier transform of the 
A*^*^) function defined in Eq. (|2.4p . and we shall see later 
that these factors ensure that the (mean-field) critical 
point occurs at 77^ = 1. The normalization factor Zi is 

defined to be (^^w)^^^x(^o)/M!. The calculation 

for Zi is straightforward, and is given in Appendix El 

The Deam-Edwards distribution can be understood as 
arising from a realistic vulcanization process in which the 
links are introduced simultaneously and instantaneously 
into the liquid state in equilibrium. Specifically, it incor- 
porates the notion that all pairs of particles that happen 
(at some particular instant) to be nearby are, with a cer- 
tain probability controlled by the link density parameter 
77^, linked. Thus, the correlations of the link distribu- 
tion reflect the correlations of the unlinked liquid, and it 
follows that realizations of links only acquire an appre- 
ciable statistical weight if they are compatible with some 
reasonably probable configuration of the unlinked liquid. 



The factor 



,M 



/Ml in the Deam-Edwards dis- 
tribution introduces a Poissonian character to the total 
number M of links. These links are envisioned to be the 
product of a Poisson chemical linking process. The fac- 
tor Z^iVo) assures that the probability of having a given 
random realization of links is proportional to the statis- 
tical weight for, in the unlinked liquid state, finding the 
to-be-linked pairs to be co-located in the liquid state to 



within the shape function exp I 



As a result of the Deam-Edwards distribution, the 
mean number of links per particle is given by [M]/N = 
rf /2. Thus, = 2[M]/N is the mean coordination num- 
ber, i.e., the average number of particles to which a cer- 
tain particle is connected, the factor of 2 results from the 
fact that each link is shared by two particles. For a de- 



tailed discussion of the Deam-Edwards distribution, see 
Ref. [lEIil. 

By using this distribution of the quenched disorder, we 
can perform the disorder average of the Helmholtz free 
energy via the replica technique, thus obtaining 

[F] ^ Y.T'{x)F^{V) 

X 

= -kBTY,T^ix)l^Z^{V) 

X 

= -knT lim Vp(x) ^"^^^"'V (2.8) 

n-!-0 ^ — ' n 
X 

We now insert the Deam-Edwards distribution to get 



[F] 



-ksT lim > 



M 



Zx(y^) 



Zxivy 



M\Zi 



(2.9) 



This disorder-averaged free energy differs from the form 
traditionally obtained via the replica technique, in that 
there is an extra replica Z^(Vo), which originates in the 
Deam-Edwards distribution. We shall call this extra 
replica the 0*^ replica, and note that it represents the 
preparation state of the system, [s^ 

The summation over the realizations of the quenched 
disorder x can be performed, following the calculation in 
Appendix \^ thus we arrive at the form 



[F] 



-ksT lim — f 



-s-o n V Zi 
which can also be expressed as 



d 

[F] = -ksT lim — InZi+n , 
n->o an 



(2.10) 



(2.11) 



where 



^ E til ^x(^o)^x(^r 



oT 



^.(^o)Z.(n'(exp(-^EnA^°'(K-s"l))) • (2-12) 



Ha 



^^^^0 li^j "=0 / 1+n 
I 



Notice that, here, the preparation state (i.e., 0*^ replica) taneously in a liquid state of fixed volume and thus the 
has a fixed volume Vq because, for convenience, we have pressure is fluctuating, whereas the measurement states 
assumed that the linking process was undertaken instan- (replicas 1 through n) are put in a fixed-pressure p envi- 
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ronment, the volume V of which is allowed to fluctuate. 
In the latter parts of the paper we shall set the pressure 
p to be the average pressure measured in the preparation 
state at volume Vq. In particular, for a given volume of 
the liquid state in which the links are made, the average 
pressure is given by 



P 



dFciVo) 



dVo 



(2.13) 



where we have introduced the Helmholtz free energy of 
the unlinked liquid FciVo) = -keTln Zc{Vo). We sup- 
pose that the excluded- volume interactions are so strong 
that the density fluctuations are suppressed, and the den- 
sity of the unlinked liquid is just N/Vq. [40] Thus, the 
mean-field value of Helmholtz free energy in the unlinked 
liquid state is 



FciVo) 



-NkeT \nVo 



2Vn ■ 



(2.14) 



Therefore, the mean pressure in the unlinked liquid state 
is given by 



P = 



Vo 



2Vo' 



(2.15) 



from which we can identify — by the standard way in 
which the second virial coefficient B2 appears in the free 
energy of the equation of state for a fluid — that B2 for 
the unlinked liquid is j^kBT ■ Thus, without having 
actually performed a cluster expansion, we see that the 
Dirac delta-function interaction with coefficient in- 
deed leads to a virial expansion with a suitable excluded 
volume, viz., v"^ / [ksT). As mentioned above, we shall 
apply this pressure in the measurement states (described 
by 1"^ through n"^ replicas), and let their volumes V 
fluctuate, in order to obtain an elastic free energy that 
can describe volume variations. In particular, by choos- 
ing the pressure p to be exactly the mean pressure of the 
liquid state, we shall obtain an elastic free energy that 
takes the state right after linking^ which has the same vol- 
ume Vb as the liquid state, as the elastic reference state. 
This issue of the state right after linking and the elastic 
reference state will be discussed in detail in Section llll Al 

In light of this construction of the pressure ensemble, 
we have the capability of learning about the bulk mod- 
ulus of the system, and to characterize volume changes 
caused by linking, a process that has the effect of elimi- 
nating translational degrees of freedoms. 

To establish an appropriate statistical mechanics for 
the fixed-pressure ensemble, we shall make the following 
Legcndrc transformation of the Helmholtz free energy, 
which leads to the Gibbs free energy G{p, T): 



P 

Q{p,T) 



dF{V,T). 



(2.16a) 
(2.16b) 



dV 1^' 
F{V,T)+pV. 

In Eq. (|2.16b[) the volume V takes the value (in terms of 
p) that satisfies Eq. (|2.16ap . i.e., the volume that mini- 
mizes the Gibbs free energy at a given pressure p. 



In the following sections, we shall first calculate 
the disorder-average of the Helmholtz free energy, and 
then make this Legendre transformation to obtain the 
disorder-averaged Gibbs free energy. This will allow us 
to explore the elasticity of the RLPM in detail. 



C. Field-theoretic description of the RLPM 

We shall use field-theoretic methods to analyze the 
disorder-averaged free energy [F] and, more specifically, 
the replicated partition function Zi+„. To do this, we 
introduce a joint probability distribution for the particle 
density in the replicated space, i.e., the replicated density 
function 

N n 



where x = [x'^ ,x^, . . . ,a;") is a short-hand for the (1-f 71)- 
replicated position c?-vector. For convenience, we intro- 
duce a complete orthonormal basis set in replica space 
{e"}S=oi terms of which a vector x can be expressed 
as 



a=0 



(2.18) 



Note that the components a;" are themselves d-vectors. 
With this notation, the density function of a single replica 
a is given by Qp^e'^ j which is the Fourier transform of the 
Q{x) field, Qp, with momentum nonzero only in replica 
a, corresponding to integrating over the normalized den- 
sities in other replicas in real space. 

The replicated partition function (|2.12|) can be written 
as a functional of the replicated density function Q in 
momentum space as 



Zi 



+n 



N . n N 

^0 i=i 



HQlQp] 



-, (2.19) 



^ a=Oj=l 



with 



HqIQp 



p 



2V 



Qpe^Q-peO 

P 

n 



p a — 1 



where the factor 



P 



,(0)\l+»,-a^|p|V2 



(2.20) 



(2.21) 



is the replicated version of the Fourier transform of the 
function A(°)(x), defined in Eq. (gH). The summation 
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2Vo 
2V 



2Fn 



2V 



denotes a summation over all momentum d-vectors with the renormalizcd coefficients 
p", one for each replica, with p taking the values p = 
X]a=o • The cartesian components of the p" take the 
values 21:171/ where L is the linear size of the system 
and m is any integer. Similarly, summations over d- 
vectors p include components having the values 21:711/ L. 

The first term on the right hand side of Eq. (j2.20p 
arises from the attractive, link-originating, interaction 
part [see Eq. (|2.12p ]: the next two terms represent the 
excluded- volume interaction in iJo, for the 0*^ replica 
and for replicas 1 through 71, respectively. 

The excluded-volume interaction is taken to be very 
strong, and thus the density fluctuations in any single 
replica are heavily suppressed. This means that Qpc-e" 
are very small for all 7^ 0, and Qpogojpa^Q = 1, cor- 
responding to a nearly homogeneous particle density. To 
manage this issue, we separate the replicated space into 
a Lower Replica Sector (LRS), in which the density fluc- 
tuations are suppressed, and a Higher Replica Sector 
(HRS), which captures the correlations between differ- 
ent replicas, and develops an instability at the liquid- 
to-random-solid transition. The definitions of the LRS 
and HRS arc (in momentum space) as follows: if two 
or more components of a replicated momentum vector 



2F"A[|°' 



(l+n) 



2T/"A 



(0) 



(2.23) 



We suppose that 



keTV 



^ rf (i.e., the excluded-volume 



repulsion is very strong, relative to the attractive effects 
of the links) , so these coefficients " are always pos- 

itive and large, relative to the energy-scale of the HRS 
that we are interested in. 



p = (p'^,p^, . . . ,p") are non zero then p is an element of 
the HRS; on the other hand, if p has zero or only one 
component p" being nonzero, and all other p^ = 0, then 
p is an element of the LRS. In addition, the LRS can be 
separated into a IRS part, in which vectors p has exactly 
one component being nonzero, and a ORS, which con- 
sists of only p ~ 0. With this separation we can rewrite 
the effective Hamiltonian as 



The interactions in Eq. (|2.22p can be decoupled using 
a Hubbard-Stratonovich (HS) transformation (for details 
see Appendix [B]). Thus, we arrive at a field-theoretic 
formulation of the replicated partition function, in terms 
of the order parameter field il: 



Hq[Qp] — - 



N7fkBT 



2F"A 



E QpQ-pAj 



(l+n) 



(0) 

peHR3 



vn, 



(2.24) 



N' 



t'^(p) E] Qpe"Q-pe° 1 (2.22) where the effective Hamiltonian is given by 



a=l 



and 



Hn[np,npec^] = ^^^^(^) X! ^p^-p^^l^'"'^ + ^^^vl{p)VLp^on_p^o 



2Vo 



^0 peHRS 



(2.25) 



dc° TTdc^exp 

Vo Jv^^o 



T/nA(0) ^ ^ P 



iN 



Ei>o(p)rjpeO e 



n 



VkeT^ 



p 



a=l 



(2.26) 
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The form of this HS transformation [see Appendix [B] especially Eq. (|B5|) ] ensures that the mean value of the order 
parameter field fl is related to the mean value of the replicated density function field Q as 

HRS: {Qp)hq ^ {np)Hn, (2.27a) 

LRS: l{Qpe-)HQ ^ {np,.)Hn, (2.27b) 

where the averages on either sides are defined via 

1 /• ^ „ n N 

^« ^ ^ / Il'^c" / n II'^S" e-^^---, (2.28a) 

:---)i/n = ^ / ^^^^pII^^p ^ ^1?^.... (2.28b) 



Q = 



The leading-order terms in iJn[f^p, i^pe"] can be constructed by expanding the InSo term in Eq. (|2.25p in powers of 
the fields flp and flp^^a , and thus we can obtain the leading-order terms in the Landau- Wilson effective Hamiltonian. 
To leading order this expansion gives 

+ ^ ^5]f^pe=f^^pe=(^l + -^^j +0((f7^)3,(17^,„)3). (2.29) 



For the LRS fields ilpe" we see the coefficients of the 
corresponding quadratic term are always positive (given 
that i/q, > 0, i.e., the excluded-volume repulsion is 
very strong), so this sector of the field theory does not 
undergo an instability. Furthermore, because these co- 
efficients (the masses, in particle-physics language) are 
very large [see Eq. (j2.23p ]. the fluctuations of these LRS 
fields Hp 

ignore these fluctuations and, for all a 
take 



The Landau theory of the vulcanization transition [26[ 
can be recovered by making the expansion of this HRS 
Hamiltonian that keeps only the leading-order terms in 
the order-parameter and the momentum p. Up to an 
additive constant and an appropriate rescaling of the or- 
der parameter, this expansion reads 



are heavily suppressed. For this reason, we Hnina] = - V (r + jpl^^r^fif] ^ 



0, 1, . . . , n, we 



pe" Ip— 



ripe" IpT^O — 0, 



(2.30) 



as a hard constraint. 

Having implemented this constraint, we arrive at the 
HRS Hamiltonian (the full form, not just the leading- 
order expansion): 



2V^"A 



(0) 



(l+n) 



where 



p£HRS 

iVfcsTlnSo, 



dc" TT<ic"cxp 



(2.31) 



(0) 



peHBS 



(l+")gjp-c 



(2.32) 



p£HRS 



^ ^ Q,p-^Q,p2^-pj^—P2 , (2.33) 



puP2£HRS 



where the potential of the links A^^^"'' has been 
momentum-expanded. This is precisely the form of the 
Landau free energy that was constructed via symmetry 
arguments in Ref. [26| . In the limit n — 0, the coeffi- 
cients become 



r oc 77^(1 — ?7^), 
V oc (??^)^. 



(2.34) 



It is straightforward to see that the r term leads to 
an instability for the link density parameter 77^ larger 
than the critical value 77^ = 1, and the lowest unstable 
modes are long- wavelength modes (i.e., p — >■ 0). [One 
should, however, keep in mind that the component p ~ 
itself, which is the ORS, is excluded from this HRS-only 
field theory; see Eq. (|2.30p .] This instability corresponds 
to the liquid-to-soft-random-solid transition, because the 
liquid state corresponds to the Hp = (in the HRS) state. 
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and becomes unstable when the link density parameter 
exceeds 1. 

There are two points that we would like to discuss 
about this Hamiltonian. First, it can be seen from the 
expansion in Eq. (|2.33p that near the critical point the ex- 
act form of interaction is irrelevant, because we have only 
kept terms to |pp in aI^^"-*, and this governs the long- 
distance physics. We have used the Gaussian-chain po- 
tential (|2.2p in this calculation. It is clear that if we were 
to change to a different potential, such as a finite-rest- 
length spring ULdrl) = f (|?'|— ^) , the long-distance (i.e., 
small-momentum) physics would be unchanged. Second, 
as we shall see in Sec. |Vl neither do the statistics of the 
elastic modulus depend on the details of the interaction 
potential; instead, they only depend on the density of 
links. This results from the fact that the elasticity orig- 
inates in the entropy of the network. The critical point 
for the vulcanization transition, and thus the entropic 
rigidity in the presence of thermal fluctuations, occurs at 
the same point as connectivity percolation does, rather 
than at the rigidity percolation critical point. The fact 
that we have a shear modulus scaling as T results from 
the entropy of the network, and not from the factor of 
T in the Gaussian-chain potential, is clear because if we 
were to change to another attractive potential, the same 
results would hold for the long-distance physics. 



D. Mean-field theory of the RLPM 

To understand the physics of the order parameter in 
vulcanization theory, and thus obtain the form of the 
stationary value of the order parameter, we need to recall 
the property of the HS transformation, Eq. (|2.27ap , which 
relates the average of the order parameter field Q, to the 
average of the replicated density function Q. According 
to this relation, we have 



1 ^ 



5\ 



{(l+n)d) 



p.O 



(2.35) 



Here, the removes the ORS part of rip. Equiva- 

lently, in real space we have 

im)Ha={^f2 '5"'+"^'^ (i - c,)) - ^,(2.36) 



which can be interpreted as 



HQ VoV" 



{n{x)) 



Hn 



1 " 

l5:[(<5('^)(x°-c°))(j('^)(zi-c]))... 



x(<5W(x"-c7))] (2.37) 



(denoted by (■•■)) with a common given realization of 
disorder one then forms the product over all replicas, 
and finally one averages over all realizations of disorder 
(an average denoted by [•••])• This interpretation can 
be understood from the definition of Hq via Zi^n, as in 
Eq. (|2J9l) . RecaU that Zi+„, as defined in Eq. ([2T2|) . 
contains thermal averages of the (1 -|- n) replicas, rep- 
resented by the factor Z^(Vb) Zy,(y)^, together with an 
overall disorder average. This validates the interpreta- 
tion given in Eq. (|2.37p . For a strict proof, see Ref. 0] . 

The structure of Eq. (|2.37p allows us to relate the 
value of the order parameter f2 to measurements on the 
system. In the liquid state, the single-particle densities 
(^g{d)^^a _ ^a-j^ each rephca are simply 1/V (or I/Vq 

for the 0*^ replica), and thus the order parameter fl van- 
ishes. In the soft random solid state, it is hypothesized 
that a finite fraction Q of the particles become localized 
around random positions. This happens when the den- 
sity of links exceeds the percolation threshold, and the 
particles that constitute the infinite, percolating cluster 
become localized. 

In the language of replica field theory, a localized 
particle remains near the same spatial position in each 
replica[41[. This is because, as we discussed earlier, each 
replica corresponds to a copy of the same disordered sys- 
tem but with independent thermal fluctuations, and for a 
particle localized as a part of the percolating cluster, it 
fluctuates around its fixed mean position, which is com- 
mon to all replicas. According to these considerations, 
it is reasonable to hypothesize the following form for the 
stationary value of the order parameter in real space: 



r2sp(x)= 



drV 



xe 



dz 



Q 



VoV- 



,(2.38) 



To arrive at this form we have assumed that the frac- 
tion of particles that become localized is Q, and the lo- 
calized single-particle density function is proportional to 
g--2 |a: -c^l replica a{= 1, . . . , n), where (^z is the ran- 
dom position near to which the particle is localized, and C 
corresponds to the uniform contraction of the entire vol- 
ume in the measurement state with respect to the prepa- 
ration state, due to linking [see the discussion following 
Eq. (|2.5ip ]. Correspondingly, this particle was near the 
position z in the preparation state, so, for replica the 
corresponding single-particle density function is propor- 
tional to e"^!^ . For convenience of notation we de- 
fine z = {z, (z, (z, . . .) as the mean position vector in 
the replicated space. The contraction ( is related to the 
change of volume as 



V_ 
Vo 



(2.39) 



This average consists of the following two steps: One first 
constructs independent thermal averages in each replica 



The localization of the particle is characterized by the lo- 
calization length ^, although for notational convenience 
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we exchange this variable for the inverse square localiza- 
tion length T = 1/^^. Because the network is heteroge- 
neous, the particles can have widely different localization 
lengths. This heterogeneity is characterized by the dis- 
tribution V{t). 

We can also write this stationarity-point order param- 
eter in momentum space: 



xe-^-f "--^ f ° ■(C-)-q45+")'') (2.40) 

The parameters that characterize this order parameter, Q 
and 'P(t), have been obtained by solving the stationarity 
condition for the Hamiltonian: 



SHn 



= 0. 



(2.41) 



The form of the order parameter (r2sp)j5 given in 
Eq. (|2.40p exactly solves the above stationarity condi- 
tion, thus we arrive at self-consistency equations of Q 
and 'P(r). In particular, the equation for Q is 



For all values of ry^, Eq. (|2.42p has a solution Q = 0, 
corresponding to the liquid state. However, for 77^ > 1, 
an additional root appears, emerging continuously from 
Q = at 77^ = 1, and describing the equilibrium amor- 
phous solid state. In Fig. [6] we show the dependence 
of the localized fraction Q on the link density, which we 
characterize by 77^ . The critical point r]^ — 1 corresponds 
to mean coordination number z of 1, and this agrees with 
the classic work on the statistical properties of random 
graphs by Erdos and Renyi [s^. For a detailed discus- 
sion and for the stationary-point distribution of inverse 
square localization lengths 7'(r), see Refs. [sl. [25|. 



(2.42) 



The contraction (, which is relevant to the elasticity of 
the random solid state, can be investigated by inserting 
the form (|2.40|) of the order parameter into the Hamilto- 
nian Hq, Eq. (|2.3ip . which yields the dependence of the 
Hamiltonian on the parameters Q, V{t) and ^. Through 
a lengthy derivation, and by keeping terms to 0{n), we 
arrive at the following Hamiltonian for the stationary 
point (cf. Appendix \C\\ : 



(SP) ^ i>g(0)iV^ ni?^{Q)N^ 



2V 



-~ NkgT \nVo - nNksTlnV + nNkBTi 9 



-(ln(2^) + C')-lny 



ln(l + l 



E 



In 

ri , . . . , r-m '^1 



Tl • • • T„ 



(2.43) 



Here, the variable f is defined as f = + a-^) , where 
T is the inverse square localization length t = 1 /^^ , and 
we have introduced the shorthand for the integrals = 
/ drP(T). The dimensionless factor B in Eq. (|2.43|) is 
given by 

= -^ + 772Q-l + e-"'Q. (2.44) 

We shall see in Section |V] that Q also controls the mean 
value of the shear modulus, as well as the amplitude of 



the disorder correlators that involve the residual stress 
fields. 

To obtain the disorder-averaged free energy, we shall 
make the stationary-point approximation: 

jj(SP) 

Zi+„~e~^. (2.45) 

Thus, we can obtain the Helmholtz free energy using 
Eq. p.lOp . arriving at the result 



[Fsp] = -fc^Thmif^" 



- 1 



-NkBT(l -^ + 0j\nV- ATfcsTy (ln(2^) + + 



„^N^ r]^NkB 



In A 



(0) 



-NkeT 



If 



2 2 



In 



T2 ksTJ 2 



'E 

m— 1 



In 



(2.46) 
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where we have used the mean- field value of Zi, from Eq. (jA3[) . and we have also made an expansion for small n of 
the renormalized excluded- volume parameter v'^ , using 



A(i+") = (A^°V+" = A^°)(l + nlnA(°^ +0(^2)). 



(2.47) 



In order to study elasticity, we shall need to know the disorder-averaged Gibbs free energy [Q], which is given by a 
Legendre transformation, Eq. (|2.16p : 



We can insert the pressure p, given by Eq. (|2.15p . and drop the slowly- varying \nV term given that ~ 77^/2 — 1 
provided the system is not close to the critical point jy^ = 1. In the limit ^ if' we arrive at 



2 AT2 r 



SPJ 



2Vo 
+NkBT 



Vo 
2 2 







ln(l + l 



2 ^ m\ 

m— 1 



In 



Tl • • ■ T„ 



(2.49) 



Using the relation p.39p . we can obtain the stationary 
values of the contraction C, that minimizes the disorder- 
averaged Gibbs free energy by solving 



0. 



In the limit 



knTV 



^, the solution is 



C~i 



OVokBT 



(2.50) 



(2.51) 



The limit -r-^rr ^ ri is the same as the limit taken 
below Eq. (|2.23p . indicating that the excluded- volume 
repulsion is much stronger than the attractive effects of 
the links. This contraction of the volume due to the 
introduction of links at a given pressure is a result of both 
the reduction of the total number of translational degrees 
of freedom, i.e., the change of the "osmotic pressure" 
and the attractive interactions induced by the links. We 
shall see later that this contraction is consistent with a 
particular phenomenological model of a disordered elastic 
medium that we shall introduce. 



E. Goldstone fluctuations in the RLPM 

1. Spontaneous symmetry breaking 

To characterize the Goldstone modes of fluctuations 
associated with the random solid state, we shall first look 
at the pattern of symmetry breaking accompanying the 
transition to this state. 

The Hamiltonian (|2.3ip for the liquid-to-soft-random- 
solid transition has the symmetry of independent trans- 
lations and rotations of each replica. The translational 
invariancc of the Hamiltonian can be readily verified by 



making the transformation 



n{x) n'{x') = n{x) = n{x' - a), (2.52) 

where a = (a*', a^, . . . , a") represents a replicated trans- 
lation. In momentum space this transformation reads 



— 5- 



p 



(2.53) 



It is easy to check that, by inserting this transformed or- 
der parameter back into the Hamiltonian ()2.3ip and mak- 
ing a change of variables, the same Hamiltonian but for 
the field fl' is recovered. Similarly, one can verify invari- 
ancc under independent rotations O = (0°, O^, . . . , O") 
with 



^6- 



(2.54) 



The order parameter in the liquid state (i.e., £7 = 0) 
has the full symmetry of the Hamiltonian, and at the 
transition to the soft random solid state it is the sym- 
metry of relative translations and rotations between dif- 
ferent replicas that is spontaneously broken. However, 
the symmetry of common translations and rotations of 
all replicas are preserved, and this reflects the important 
notion that from the macroscopic perspective the system 
remains translationally and rotationally invariant, even 
in the random solid state. This entire pattern of sym- 
metry breaking amounts to an unfamiliar but essentially 
conventional example of the Landau paradigm. 

This broken symmetry of relative translations and ro- 
tations between different replicas can be understood as 
a result of particle localization. Because a delocalized 
liquid particle can explore the whole volume via its ther- 
mal fluctuations, and in thermal equilibrium its positions 
in different replicas are uncorrelated, the liquid state is 
invariant, under separate translation and rotation of indi- 
vidual replica. On the contrary, for a localized particle. 
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its positions in the various replicas are strongly corre- 
lated, and therefore the symmetries of relative transla- 
tions and rotations are broken. 

It is straightforward to verify that the form of the 
random-solid-state order parameter, Eq. (|2.38|) . correctly 
implements this pattern of symmetry breaking. To see 
this, we can use the complete orthonormal basis in replica 
space defined in Section III C[ and define an alternative 
basis involving a "replica (almost) body-diagonal" unit 
vector 



1 



(2.55) 



Relative to e^, we may decompose a (l + n)d dimensional 
vector X into its longitudinal (A) and transverse (r) com- 
ponents: 

x = x\+Xr, xx = ix-ex)ex, Xr = x - xx .(2.56) 

Note that xx and x^- are both (1 + n)(i-dimensional vec- 
tors, but Xx has only d-dcgrccs of freedom (given by 
X ■ ix), and Xr has only n,c?-degrccs of freedom. 

By this decomposition, the vector z = (z, ^z, ^z, . . .), 
which characterizes the mean positions of a particle in 
the stationary-point state, can be written as 



= i/l + nC^ zex, 



(2.57) 



which points purely in the direction. As a result, the 
stationary order parameter, Eq. (|2.38p . can be written as 



T \ 2 



xe 2 



VoV" 



Q / drP(T) 



D 2 1^"^ I" 



27r 



2tt 



r(l + nC2) 



(2.58) 



where in the last line we have integrated out the d- 
dimensional vector z. It is evident that this value of 
order parameter does not depend on xx, which means 
that it is invariant under translations in the ix direction, 
corresponding to common translations and rotations of 
all replicas ( albeit appropriately contracted by C, in repli- 
cas 1 through n ). This stationary order parameter is 
shown schematically in Fig. [Ija) for two replicas. The 
Gaussian-like form in the Xr direction indicates a "con- 
densation" between different replicas. This is called a 
molecular bound state in Ref. |18j . 



2. Goldstone fluctuations 

With the pattern of continuous symmetry breaking 
just outlined, we can write down the form that the order 



parameter takes when it is subject to "Goldstone fluctu- 
ations" : 



xe 



■,(2.59) 



where R^{z) 



Therefore the Goldstone deforma- 



tion of the order parameter is parameterized by the 
n independent functions {i?^(z), . . . , i?"(z)}. The sta- 
tionary form of the order parameter, Eq. (|2.38p . de- 
scribes a system in which the mean positions of the 
replicas of the thermally fluctuating particles are located 
at (a:;°, x^, . . . , x") = (z, Cz, ■ • ■ , C-^^)- We shall refer to 
these positions as the centers of the thermal cloud''! By 
comparing the undeformed order parameter (j2.38|) and 
the "Goldstone-deformed" one, Eq. (|2.59p . we see that 
the Goldstone-deformed order parameter describes a sys- 
tem in which the mean positions of the replicas of the 
fluctuating particles are displaced from (z, C,z, . . . , C,z) to 
(z, i?i(z), . . . , -R"(z)). Thus, i?"(z) (a = 1, 2, . . . , n) rep- 
resent the deformed mean positions in the measurement 
replicas. 




(a) 



(b) 



FIG. 1: (a) Schematic plot of the value of the order param- 
eter (brightness) at the stationary point, for tlie illustrative 
two replicas, labeled by and . (b) Schematic plot of 
the value of the order parameter (brightness) for a Goldstone 
deformation of the stationary point for two replicas. 



We require that the deformations C,z — > R°'{z) be 
pure shear deformations. This constraint can be ex- 
pressed as det (c)i?f(z)/9((^z)j) = 1; it guarantees that 
the Goldstone fluctuation does not excite the LRS (i.e., 
each replica still has homogeneous density), which would 
be extremely energetically costly, owing to the large 
excluded-volume interaction. The ORS has already been 
removed from the theory, and one can easily check that 
it remains zero in this Goldstone-deformed order param- 
eter. The vanishing of the order parameter in the IRS 
can be verified by taking the momentum-space Goldstone 
deformation and making a change of variables: 
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q=pe° 



dTV{T)e 



_l|L_ipi?-(^) 



Q 



dR°' V 



V Vo 



dz 



dR°' 



dTV{T) 



-ip i?" (z) 



}d) 



(2.60) 



This result indicates that the deformed state has the 
same local density as the (undeformcd) stationary-point 
state. 




R iz) 



FIG. 2: Example of a Goldstone-deformed state. The system 
of replicas through n are shown. The mean positions of the 
replicas of a thermally fluctuating particle are displaced to 
(z, R^{z), . . . , R"{z)) in this Goldstone-deformed state, which 
characterizes an n-fold replicated deformation field. Here, for 
simplicity, we only show spatially homogeneous deformations, 
and it is worth noticing that the volumes of the measurement 
replicas, i.e., replicas 1 through n, axe contracted by a factor 

We ought to clarify the following point about this 
Goldstone deformation. As we have already mentioned, 
the symmetry that is broken at the transition is that of 
relative translations and rotations of the various repli- 
cas, the symmetry of common translations and rotations 
remaining intact. As a result, the Goldstone deforma- 
tion should be constructed via z-dependent translations 
of the order parameter in the Xr direction, i.e., the bro- 
ken symmetry direction. However, if we look at the 



deformation field defined by U = R — z, we find that 
U = (0, R^{z) — C^i ^^(^) ~ C-^^: ■ • ■) is in fact not in the 
broken symmetry direction ir, because it has an xx com- 
ponent, viz.. 



Ux = iU- ex) ex = 



C 



Y,U%z)ex. (2.61) 



This Ux component is actually redundant. This can be 
seen by decomposing the quadratic form \x — z — U\'^ as 



\xx 



zx-Ux? + \Xr-Ur?. (2.62) 



and noting that in the form of the order parameter ()2.59|) 
one can change the integration variable (which is a d- 
dimcnsional vector) from z to 



1 



U-ex, 



(2.63) 



so the longitudinal component of the (l+n)(i-dimcnsional 
vector y = (y, Cy, . . . , C^y) is yx ~ zx + Ux- The Jacobian 
of this change of variables is unity, provided that each 
deformation z — > U"{z) is a pure shear deformation [4^ . 
With this change of variables the Goldstone-deformed 
order parameter attains the form 



e 2 



e 2 



Q 



VnV" 



(2.64) 



where the transformed deformation field U' is defined 
via U'riy) — lJr{z). With this change of variables, the 
order parameter field, 5lGF(i) of a;, is transformed to 
a new field, U,'Q^{y) of y, which can be viewed as be- 
ing the stationary point ^sp{y) but locally translated 
purely in the yr directions, because the deformation is 
R'{y) = y + U!^{y). Comparing with the deformation 
before the change of variables, R{z) = z + U{z), it is evi- 



dent that the Ux component is actually removed, and the 
deformation field U only affects the direction. There- 
fore, Ux is a redundant component in this field-theoretic 
description of the Goldstone fluctuation. Note that in 
these two representations of the Goldstone fluctuation, 
Eqs. (|2.59p and (|2.64p . the number of degrees of freedom 
of the deformation field, U (z) or U!^{y), is nd, because in 
U (z) one has the constraint = 0. 
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The reason we choose to adopt the form of Goldstone 
fluctuation given in Eq. (|2.59p is that in the true physical 
system that we are intending to describe, the preparation 
state (replica 0) is not deformed. Although in the field 
theory the 1 + n replicas feature symmetrically (apart 
from the contraction ^) , physically, one should only have 
Goldstone fluctuations that deform replicas 1 through n, 
as these are the replicas associated with the measurement 
states, on which deformations are actually performed. 
Therefore, although Eqs. (|2.59p and (|2.64|) are mathe- 
matically equivalent, Eq. (j2.59p provides a better physi- 
cal description of the Goldstone fluctuations j43j . 



3. Energetics of Goldstone deformations 

To obtain the energy of a Goldstone deformed state, 
we take the momentum-space version of the Goldstone- 
deformed order parameter, Eq. (|2.59p . 



xe-^-*^-^(-)-Q,S(i+"''', (2.65) 



and insert it into the Hamiltonian (|2.3ip . After a lengthy 
calculation (see Appendix |D]), similar to the one for the 
stationary-point free energy, we arrive at the energy of 
the Goldstone deformed state: 



r(SP) 



(2.66) 



Here, we use the short-hand 

■^r{zi,Z2) = {R{zi)-R{z2)f-{l + n){zi-Z2f 

n 

=J2{{R-{z,yR-{z2)r-{z,~Z2)^ (2.67) 



(SP) 

to denote the deformation. The term H^^ is the Hamil- 
tonian at the stationary point, as given in Eq. (|2.43p . and 
the term i?*^ accounts for the increase in the energy due 
to Goldstone deformation, and is given by 

^ j dzidz2Ki{zi, Z2)^ r{zi, Z2) 

-i— [ dzidz2dz3dz4K2{zi,Z2,Z3,Z4) 
^bT J 



x^'^(zi,Z2)*fl(z3, ^4) -riArfcBTy(C'- 152.68) 

where we have kept terms to linear order in n and 
quadratic order in V^lj in this expansion. The functions 
Ki(zi, Z2) and K2{zi,Z2, z^,Z4) are bell-shaped functions 
of the distances (zi — Z2) and (zi — Z2), (^2 — 23), 
(2:3 — Z4), as shown in Fig. [S] They are independent 
of the centers of mass, (zi + Z2)/2 for Ki{zi,Z2), and 
(zi -I- Z2 -I- Z3 -I- Z4)/4 for K2{zi,Z2, Z3, Z4). The forms of 
Ki{zi, Z2) and K2{zi, Z2, Z3, Z4) are given in Appendix ID] 



In specific, _fi'i(zi,Z2) has the following schematic func- 
tional form 



Ki{zi,Z2) 



drC(r 



(2.69) 



which is a superposition of Gaussian distributions on the 
scale of r, where r is a certain combination of localization 
lengths ^, weighted by the distribution of those localiza- 
tion lengths. The kernel i<'2(-Zi, 2:2, Z3, Z4) has a similar 
structure, except that it also contains delta-function fac- 
tors such as 5{zi — Z3), which should be associated with 
the scale of the short-distance cutoff of the theory, which 
is also on the scale of typical localization length. 

7, V. 

Kl{zi,Z2) = Z 



A'2 (21,2:2,23. 24) = z 




+ z 



FIG. 3: Diagrams for the functions Ki{zi,Z2) and 
K2{zi, Z2, Z3, Z4). In these diagrams, the straight lines rep- 
resent delta functions, and the wavy lines represent Gaussian 
potentials between pairs of points, averaged over the distri- 
bution of localization lengths, as indicated schematically by 
the summation signs, which also indicate symmetrization over 
the arguments. The typical localization length provides the 
characteristic lengthscale for Ki and A'2. The full expression 
is listed in Appendix iDl 

The form of the energy increase due to Goldstone fluc- 
tuations, Eq. (|2.68p . can be understood intuitively as fol- 
lows. This energy actually describes the elastic energy of 
replicated shear deformations of the system, as we have 
explained in Section [HE 21 (i.e., the Goldstone fluctua- 
tions in the RLPM are replicated shear deformations). 
Thus, it is evident that the first term in Eq. (|2.68p repre- 
sents a coupling of the deformations i?(zi) and i?(z2) at 
the points zi and Z2. The coupling function, Ki{zi, Z2), 
given in Eq. (|D10p . has a magnitude controlled by the 
probability Q for a particle to be localized, and a length- 
scale controlled by the typical localization length; these 
two quantities are, in turn, determined by the link den- 
sity parameter 77^, which is the control parameter for the 
RLPM. In particular, the first term in Ki{zi,Z2) car- 
ries a factor Q^, which is the probability for both of the 
two mass-points, zi and Z2, to be in the infinite clus- 
ter, and therefore to have an elastic interaction between 
their deformations. The second term in Ki(zi, Z2), which 
involves a summation over m from 2 to 00, takes into ac- 
count the interactions between these two points that are 
mediated via other mass points. The four-point term. 
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with couphng function K2{zi, Z2, Z3, Z4), is a bit more 
compUcated. It couples repUcas of the shear deformation 
to one another and to themselves, also on the scale of the 
localization length. These two terms in the elastic energy 
embody the statistics of the elastic free energy of the soft 
random solid in the language of replicas, and thus encode 
information about the statistics of the quenched random 
elastic parameters. In Section IIIII we shall uncover this 
statistical content via the introduction of a phenomeno- 
logical model inspired by this elastic energy. 

In order to compare with the phenomenological elastic 
free energy that we shall discuss in Section lllll it is useful 
to make an alternative decomposition of the Hamiltonian 
to the one made in Eq. p.66p . That decomposition ()2.66|) 
was in terms of the stationary-point part and the part due 
to fluctuations. The alternative decomposition reads 



{R) 



(2.70) 



The relation between the two decompositions is given by 



(0) 



{R) 



(SP) 

n 

*i7 



HO, 
-HO, 



(2.71) 



where HO accounts for the energy of the stationary 
point, measured with respect to the state right after link- 
ing, which is actually the elastic energy of the contraction 
C, i.e.. 



HO 



NkBT^iO'-^)- (2.72) 



In Eq. p.70p . we are separating the Hamiltonian into 
two parts. H^'' gives energy of "the state right after 
linking" which is a state that has not been allowed to 
contract after the links were made and thus has the 
same volume and shape as the liquid state. In the 
state right after linking, the mean positions of the repli- 
cas of the thermally fluctuating particle (i.e., the cen- 
ters of the thermal cloud) are located at the positions 
(a;°, cc^, . . . , a;") = {z,z,...,z). The other part of the 
Hamiltonian is the elastic energy of the deformation away 
from this state, . This is different from the separa- 
tion made in Eq. (|2.66p . in which one has the stationary 

(SP) 

point energy and the energy-increase due to Gold- 

stone deformations i?*". 



III. PHENOMENOLOGICAL APPROACH TO 
THE ELASTICITY OF SOFT RANDOM SOLIDS 

A. Phenomenological nonlocal elastic free energy 

As discussed in Section HI in the classical theory of 
rubber elasticity [Ij, rubbery materials are modeled as 
incompressible networks of entropic Gaussian chains, and 
the resulting elastic free energy density is given by 



for spatially uniform deformations r — ;> A-r. Incompress- 
ibility is incorporated via the constraint detA = 1. For 
the shear modulus fi, the classical theory gives the result 
?ic ksT, where ric is the density of effective chains in the 
network. 

The phenomenological model that we now discuss is 
in the same spirit as the classical theory of rubber elas- 
ticity. However, to account for the heterogeneity of the 
medium we need to introduce the additional feature of 
quenched randomness into the model, and thus the en- 
tropic Gaussian chains are allowed to be of heterogeneous 
length and density. Furthermore, the classical theory is 
a local elasticity theory, which is valid at length scales 
that are much longer than the effective chain length of 
the polymers. By contrast, our phenomenological model 
is a nonlocal theory, which explicitly takes into account 
the finite length of polymer chains, as well as their vari- 
ations. 

Inspired by the form of the energy of Goldstone fluctu- 
ations determined from the RLPM in Section IIIE3[ we 
choose the following elastic free energy Fq, associated 
with a deformation of the soft random solid state that 
maps the mass point at z to the new location R{z): 

Fg = i I dz, dz2 G(zi,Z2)(|i?(zi) - i?(22)|'- \Z1 - ^2P) 



Ao 

' 2 



(3.2) 



/ 



■ Tr A ^ A 



(3.1) 



where G{zi,Z2) is a nonlocal harmonic attraction that 
serves to pull the two "mass points" (i.e., coarse-grained 
volume-elements) at z\ and Z2 towards one another. The 
kernel G{zi,Z2) originates in the entropy of the molecu- 
lar chains of the heterogeneous network, and we model it 
as "zero rest-length" springs having random spring coef- 
ficient. Notice that 0(^1,22) is a coarse grained conse- 
quence of many molecular chains and, more importantly, 
is an entropic effect and does not depend on the choice 
of precise form of microscopic attractive interactions. 

We take G(zi, Z2) to be a quenched random hmction of 
the two positions, zi and Z2, symmetric under zi ^ Z2- 
We assume that the disorder average of G{zi,Z2) is 
G'-^^zi — Z2) = [G{zi, Z2)], i.e., is translationally in- 
variant. Furthermore, we define the fluctuation part of 
G(zi,Z2) to be G(i)(zi,02) = G(zi,z2) - GW(zi - Z2). 
In the following analysis, we assume that G(i) < G(o) m 
order to make a necessary perturbative expansion. 

In the second term in Eq. (|3.2p . the determinant of the 
deformation gradient tensor A.y (z)[= dRi/dzj] captures 
the change of the volume and, correspondingly, the pa- 
rameter Ao, which we take to be large, heavily penalizes 
density variations. This large Aq results from a competi- 
tion between (i) repulsions (either direct or mediated via 
a solvent, e.g., excluded- volume) , and (ii) intermolecular 
attractions and external pressure. 

In discussion of elasticity that follows, we exploit the 
notions of a "reference space" and a "target space" for 
any deformation R{z). The reference space, labeled by 
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the d-dimensional vector z, is the space before the de- 
formation, whereas the target space, labeled by the d- 
dimensional vector R{z), is the space after the deforma- 
tion. 



B. Disorder average of the phenomenological 
model via the replica method 

To make a comparison with the RLPM, and thus to 
obtain information about the statistics of the nonlocal 
kernel G that characterizes the disorder present in the 
phenomenological model, we shall use the replica method 
to average the elastic free energy p.2|) over the quenched 
disorder, whose statistics will be specified below. 

Wc follow a recipe similar to the one used in Sec- 
tion |TlB] [sec Eq. dill)]. The elastic free energy, Eq. 
contains the random ingredient G. As with the RLPM, 
the physical quantity to be disorder-averaged is the free 
energy at a given pressure, but now with the deforma- 
tions R{z) as thermally fluctuating field. Therefore, we 
need to take T, defined in Eq. p.2p . as the effective 
Hamiltonian, because it is the elastic energy for a given 
deformation field specified by R{z), and then calculate 
the free energy at a given temperature, via the partition 
function 



where we have used [...]= J VG Vg ■ ■ ■ once again to 
denote a disorder average, but this time over the values 
of the quenched random kernel G, weighted by an as-yet 
unknown distribution Vg- In the present setting, we do 
not have a "zerotli replica" as such a replica arises from 
the Deam-Edwards distribution of the links, and this is 
not the type of quenched disorder that we have in mind. 
Rather, in the present setting we regard the distribution 
of disorder Vg as a physical quantity that is unknown 
but is to be determined through a comparison with the 
analysis of the Goldstone fluctuations of the RLPM. The 
replica partition function is then given by 



[Zg] = j VGVgZ^ 
JvGVgJ flVR^e-^'^ 



Jl VR 



a = l 



^ 1 



(3.5) 



Zg = / Pi?e-ro(fl(-))A-i.T 



(3.3) 



with r depending on the quenched randomness through 
its kernel G. The Gibbs free energy is related to this par- 
tition function via Q = —ksTlnZ, and Q is the quantity 
that should be averaged over the quenched disorder. Note 
that it is the Gibbs free energy, instead of the Helmholtz 
free energy, that is related to this partition function Z , 
because in the elastic energy one has a fixed pressure, 
which is accounted for by the Aq term, the volume being 
allowed to fluctuated. The disorder average of the Gibbs 
free energy can be computed using the replica technique: 



in which we functionally integrate over the configurations 
of the n-fold replicated displacement fields R". We have 
also introduced the effective pure Hamiltonian r„ that 
governs the replicated deformation fields: 



r„ = -A:BTln[( 



-ES=irG(fl°(^))/feBT] 



(3.6) 



-keT / VGPG\nZG 



-keT I VGVg lim 



'ksT lim 



dn 



ri-i-O n 



\Zg 



(3.4) 



The exponential and the logarithm in Eq. p.6|) can 
jointly be expanded in terms of cumulants, and thus we 
arrive at the form 



n n 

kBT{-[Y,^G{R"{z))/kBT]^ + -[ J2 rG{R"{z})rG{R^{z))/kBT]^----} 



a, 13=1 



(3.7) 



where [■ . ■]c are connected statistical moments (i.e., cumulants) associated with the probability distribution of the 
disorder Vg-, and the omitted terms are 0[{T /ksTY]. The elastic energy Tq for a given realization of disorder G and 
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a given deformation field R{z) is given in Eq. p.2p : inserting tliis form for Tq we obtain 

"'^ a=l "'21,22 

[G{Z1, Z2)G{Z3, Z4)]c'ifRizi, Z2)'^B.iz3, Z4) + ©(^fj), 
I 



1 



(3.8) 



and we remind the reader of the definition of Vt'jj, first 
given in Section lilE3l 



■^R = (i?(zi) - R{Z2)Y - (1 + n)(zi - Z2f 



(3.9) 



Up to quadratic order in "^fji, the effective pure Hamil- 
tonian r„ of Eq. p.Sp has precisely the same form as 
the energy of the Goldstone fluctuations (|2.68p . derived 
microscopically from the RLPM. Thus, the RLPM actu- 
ally provides a derivation of the phenomenological model 
we proposed in Section [ill Al and justifies, from a micro- 
scopic perspective, the phenomenological elastic free en- 
ergy p.2p with its quenched randomness. Therefore, the 
probability distribution Vg of the quenched randomness 
in Eq. p.2p is contained in the RLPM. By comparing 
the two schemes, i.e., Eqs. (|2.68|) and (|3.8p . we arrive at 
a statistical description of the quenched random kernel 
G in the phenomenological model (|3.2p . as we shall now 
show. 



C. Comparing the Gibbs free energies of the 
RLPM and the phenomenological model 

The RLPM is a semi-microscopic random network 
model, our analysis of which led to the disorder-averaged 
Gibbs free energy (|2.11I2.16|) : 



[S] 



d 

-ksT lim — InZi-f 
n-s-o On 



with 



(3.10) 



(3.11) 



The functional integration here is over all possible con- 
figurations of the order-parameter field Q,. By contrast, 
in the phenomenological model the replicated partition 
function Z„ involves a functional integration over the 
n-fold replicated deformation field i?", as in Eq. (|3.5p . 
To obtain the equivalence between the RLPM and the 
phenomenological model it is useful to proceed in two 
steps. First, we note that in the random solid state the 
Boltzmann weight in Eq. (|3.1ip [i.e., exp (— TJo/fcsT)] 
is heavily concentrated near the stationary point i^sp 
and the Goldstone-deformed states Oqf, and decreases 
steeply for other sectors of fluctuations. This suggests 
that we parametrize the fluctuating field f2 in terms of an 



"amplitude" which we take to be il^p, plus "radial" fluc- 
tuations around it, together with an appropriate set of 
generalized "angular" Goldstone variables, which are the 
n independent d- vector fields in as, e.g., in Eq. (|2.65p . 
We then make an exact change of functional integration 
variables, from VL to these radial and angular variables, 
and this introduces a corresponding Jacobian factor. The 
second step is to recognize that the radial fluctuations 
are "massive" (i.e., they have restoring forces, in con- 
trast with the angular fluctuations, which are "mass- 
less"). Thus, if we were to integrate these radial fluc- 
tuations we would obtain small corrections to the terms 
of the remaining, effective, angular- variable theory. We 
therefore elect to treat the radial variables at the classical 
level, which amounts to neglecting the radial fluctuations. 
Under this condition, the aforementioned Jacobian factor 
does not depend on the angular variables, and therefore 
it contributes only a constant multiplicative factor to the 
functional integral, which can be safely omitted. This 
procedure enables us to arrive at the following approxi- 
mate form for the replica partition function of the RLPM, 
Eq. (IXTTI) : 



Zi 



W VR 



(3.12) 



with i?^^^^ and given in Eqs. (|2:43)) and ((2:68)) . 

In our phenomenological model, introduced in Sec- 
tion nil Al the disorder-averaged Gibbs free energy is 
given by 



with 



-ksT lim ^ 
n->-o on 



a=l 



-r,jkBT 



(3.13) 



(3.14) 



where r„ is given in Eq. 

The Gibbs free energies — for the RLPM and for the 
phenomenological model — are supposed to be equal, up 
to an additive constant, because they both capture the 
Gibbs free energy of a soft random solid system having 
elastic deformations. It is this equality that we shall now 
exploit to characterize, via the RLPM, the distribution 
of quenched disorder Vg hi the phenomenological model. 
Actually, we can directly identify the Hamiltonian r„ 
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with Hq , because the functional integration over the 
repHcated deformation field R" is common to both the 
RLPM and the phenomenological model, in the sense 
that the deformation z — > R°'{z), in both the RLPM 
and the phenomenological model, takes the state right 
after linking as the reference state. Therefore, we have 
the relation 

r„ = 4^). (3.15) 

Notice that, here, the RLPM Hamiltonian is h[^\ not 

H^"^, because it is Hq^^ that is the energy measured 
from the state right after linking, which matches the def- 
inition of reference state in the phenomenological theory, 
whereas H^'^ is the energy measured from the stationary 
point, which differs from the state right after linking by 
the energy associated with the contraction given in 

Eq. srm . 

By the comparison stated Eq. (|3.15p . we arrive at the 
following determination of the quenched-disorder charac- 
teristics of the phenomenological model (LHS) in terms 
of the elastic properties of the RLPM (RHS): 

[G{zi,Z2)]c = Ki{zi,Z2), (3.16a) 
[G{zi,Z2)G{z3,Z4)]c = K2{zi,Z2,Z3,Z4),{3.l6h) 
Ao = i^^nl, (3.16c) 

where uq = N/Vq is the number-density of the particles 
in the preparation state. The functions Ki{zi,Z2) and 
K2{zi, Z2, Z3, Z4) are defined in Appendix iDl and have 
been discussed in Section Hi E 31 



inward pressure, as well as introducing random stresses, 
as shown in Fig. U] As a result, immediately after cross- 
linking the state is not stable, but relaxes to a new stable 
state, determined by the particular realization of ran- 
domness created by the cross-linking. In the following 
discussion, we shall use the following nomenclature: 

R{z) = z ^==^ the state right after linking, (4.2a) 
R{z) — z = (z + v{z) <:=> the relaxed state. (4.2b) 

The state right after linking, here, is the same as the one 
just defined, following Eq. (|2.70p . which has energy H^^ 
in the RLPM, because they both describe the state that 
has undergone no deformation since being linked. 




(a) (b) (c) 



FIG. 4: Schematic plot of the relaxation process under a fixed 
pressure, (a) The liquid state with no linking, (b) The state 
right after linking, cross-links are added to the system, and an 
infinite cluster is formed. This state is not stable, because of 
the inward pressure and local stresses, (c) The relaxed state. 
The system undergoes a uniform contraction and random lo- 
cal deformations that release the unbalanced stress introduced 
by cross-linking. 



IV. PHENOMENOLOGICAL MODEL AT FIXED 
DISORDER: RELAXATION, EXCITATION, AND 
DEFORMATION 

A. Relaxation to the stable state at fixed disorder 

The free energy F provides a natural description of 
the heterogeneous elasticity of soft random solids. How- 
ever, its stable state is not R{z) = z (i.e., the state 
R{z) = z does not satisfy the stationarity condition 
6r/6R{z) = 0). There are two reasons for this insta- 
bility. First, the attraction G causes a small, spatially 
uniform, contraction [the fractional volume change be- 
ing 0(1/ Ao)]. Second, the randomness of G additionally 
destabilizes this contracted state, causing the adoption 
of a randomly deformed stable state. We denote this 
relaxation as 

z^z = (z + v{z), (4.1) 

in which ^ describes the uniform contraction and v{z) 
describes the random local deformation. This relaxation 
process can be understood in the setting of the prepara- 
tion of a sample of rubber via a hypothetic instantaneous 
cross-linking: cross-linking not only drives the liquid-to- 
random-solid transition but it also generates a uniform 



By writing the relaxation as z —> z [= (^z + v{z)] we 
are making the approximations that the contraction C is 
homogeneous and that the random deformations v{z) are 
pure shear, which means that any randomness in the bulk 
deformation is ignored. This can be understood by look- 
ing at the orders of magnitude of the deformations. The 
uniform contraction is of order 0(G^''-' /Aq), and the ran- 
dom local shear deformations are of order 0{G^^^ /G'^'^^). 
The random local bulk deformations is, however, of order 
0{G^^^ / Xq), and is thus much smaller than the other two 
deformations, given the assumptions that (i) the fiuc- 
tuations of the shear modulus arc much smaller than 
the mean value (the shear modulus corresponds to G, 
as we shall see later), and (ii) the shear modulus is much 
smaller than the bulk modulus. 

With these assumptions, we can insert the form R{z) = 
z(z) into the stationarity condition, and solve for the re- 
laxed state, which is characterized by ( and v{z). As we 
have just discussed, for the contraction, only the homo- 
geneous part is admitted, so the variational equation for 
C assumes G^^'' = and thus v{z) ~ and, stationarity 
requires 
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Thus, for the present model, Eq. (|3.2[) . we have 

0=^{^ J dzidz2G(")(^i,Z2)(C'-l)kl-^2p 



4/d.(c^-i)^} 
= Vi,(dpC + Ao(C'-l)dC'^-i). 



(4.4) 



By solvmg this equation to leading order in p/Xo, we 
obtain 



C « 1 - (p/dAo), 



where 



P=\j dz2 (Zi- Z2)'G(°) (Z1-Z2). 



(4.5) 



(4.6) 



As we shall see below, p is actually the mean shear mod- 
ulus. 

The stationarity condition for the random local defor- 
mation v{z) reads 



5T 



5Va{z) 



and for the present model, Eq. 
comes 



0, (4.7) 
this condition be- 



2{CZa+Va{z)) j dz2G{z,Z2) 
-2 J dz2G{z, Z2){Cz2m + Va{z2)) 



(4.8) 



Here, the last term, XQda{diVi{z)^ , is associated with 
density variations, and arises from the variation of the 
second term in the elastic free energy (|3.2p . which is 

^ / dz[det(^dRi{z)/dzj) — l] . In the following discus- 
sion we shall call this the "bulk term" in the elastic free 
energy, and we have made the definition Aq = C^'^^^Aq; 
see Appendix [El for the expansion. 

The stationarity equation (j4.8p for v{z) can be solved 
perturbatively, by assuming that G^"^ is of zeroth order 
and that G^^^ and v{z) are of first order; see Appendix|E] 
for the explicit calculation. In momentum space, the 
result is 



■ fp ■ fp 

2Dp X'o + 2Dp' 



(4.9) 



where p is a rf-dinicnsional momentum vector. The no- 
tation fp and Dp are defined as 



fa.p 



G 



(1) 



p'=0 



G, 



(0) 



Gi"'. 



(4.10) 



Notice that fa^p is actually the random force in the state 
that is contracted but not yet equilibrated for random- 
ness. The definitions of the projection operators p'^ and 
p""" are respectively. 



(4.11) 



We use bold letters to denote d-dimensional rank-2 ten- 
sors, and add an overhead arrow [such as w] to denote 
vectors, when needed. 

In the solution (|4.9p . the second term is much smaller 
than the first term, due to the large bulk modulus Ag. In 
the incompressible limit (i.e., Aq — > 00), we have 



2Dr, 



(4.12) 



which is a purely transverse field, meaning that it satis- 
fies PiUip ~ or, cquivalently, that diVi{x) = 0, which 
is the only type of deformation that can occur in an in- 
compressible medium. 



Excitation around the relaxed state at fixed 
disorder 



In order to obtain a description of the elasticity of the 
relaxed state z ^ z, which is a stable state and thus 
relevant for experimental observations, we re-expand the 
phenomenological elastic free energy (j3.2p around the re- 
laxed state. This amounts to taking the relaxed state 
z(z) [ = C,z+v{z)\ as the new reference state, and deriving 
the elastic free energy for deformations (i.e., excitations) 
relative to this state. 

To do this, we study the free energy for the following 
elastic deformation: 



z -> i?(z) = z(z) + u(z). 



(4.13) 



where u{z) is a deformation away from the relaxed ref- 
erence state. It will be convenient to make the following 
change of the independent variables: 



z z{z) = ^z + v{z), 
which has the Jacobian determinant 



(4.14) 



d. 



Zi 



dzj 



C^(l + r^9,^;,(z)). (4.15) 



With this change of variables, the phenomenological elas- 
tic free energy is expressed as 
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r = ro + - I dzidz2J{zi) J{Z2) G{zi,Z2)[\zi+u{zi)-Z2-u{z2)\ -\Z1-Z2 



^ / d~zJiz)-'[jiz)dct 



dzj + Ui{z) 
dzi 



(4.16) 



where we have made the definitions 



G(zi,Z2) = 

u{z) = 
R{~z) ^ 



G(z(zi),z(z2)), 

u(z(z)), 

z + u{z), 



(4.17a) 
(4.17b) 
(4.17c) 



with z{z) denoting the mapping of the mass point z in 
relaxed state back to the mass point z in the state right 
after linking, i.e., the inverse of the z{z) mapping. The 
change of the free energy due to choosing a different ref- 
erence state is defined as 

Tn= i Jdzidz2GizuZ2){\h - S2\^-\Z1 - Z2\^),{4..18) 

which is a constant for any given realization of the ran- 
domness. 

In order to obtain a direct description of the elas- 
tic energy relative to the relaxed state, we expand the 
quenched random nonlocal kernel in the relaxed state, 
G(zi,Z2), defined in Eq. (|4.17ap as 



(1) 



(o)^ 



where v is the random local deformation field defined in 
Eq. (|4.ip . We then expand the elastic free energy (|4.16p 
as a power series in the small deformation u{z) away from 
the relaxed state z. The computation of this expansion 
is given in Appendix IF] 

As we shall show in Section |Vl the statistics of the 
quenched randomness present in this phenomenologi- 
cal theory can be determined via a comparison with 
the RLPM. Through this comparison, we find that, the 
lengthscale of the nonlocal kernel G is actually the typi- 
cal localization length, which is small compared to the 
lengthscales on which our theory of elasticity applies, 
because the deformations in this theory are associated 
with Goldstone fluctuations in the RLPM, which feature 
lengthscales larger than the typical localization length. 

Thus, it is reasonable to make a local expansion of the 
elastic energy (|4.16p relative to the relaxed state in terms 
of the strain tensor e. The resulting form, which we shall 
call "the local form of the elastic energy relative to the 
relaxed state" is in the form of Lagrangian elasticity. As 
will be seen in Section IV Bl the advantage of this local 
form of the elastic energy is that one can extract from 
it large-distance behavior of the disorder correlators of 
the elastic parameters, which turn out to be universal. 



The calculation for this local expansion is given in Ap- 
pendix |F1 The resulting local form of elastic energy is 



dz{Tr((T(z) • e(z)) 



-/i(z)Tre(z)^ 



(Tre(z))2}, (4.20) 



where the strain tensor relative to the relaxed state is 
defined as 



g (z) ^ -(^^ + + 

~ 2\ dzi dzj dzi dzj 



.21) 



and the heterogeneous elastic parameters, viz., the resid- 
ual stress cr, the shear modulus fx, and the bulk modulus 
A, are given in momentum space, by 



An 



G 



dqidqj 

fa,p I ~ 

\v\ 

■ip- fp 



(1) 

P-9.P 



q=0 

Ja,p ( ~ T I ~ T ~\ 
I -.12 [PiPja.p ^PjPia.pj^ 



+ (4.22a) 



\P\ 



IpI' 



AoK)5p + 2{^^-pl/o<5p} 



(4.22b) 
(4.22c) 



In the expression for cTy-^p we have kept terms only to 
leading order in the momentum p (see Appendix[F]for the 
derivation) (4^ . It is worth mentioning that, to leading 
order in the momentum p, this residual stress satisfies 
the stability condition pi (Jij^p = 0, because the reference 
state of this elastic free energy, the relaxed state, is a 
stable state. This will also be shown more directly in the 
final results in Section IV Bl 



C. NonafHne deformation at fixed disorder 

Because of the quenched disorder present in the elastic 
parameter G of our phenomenological model, Eq. p.2p , 
upon the application of external stress, the system will 
respond by adopting a strain field that is nonaffine. This 
means that the strain tensor will be spatially inhomoge- 
neous even though the applied stress is homogeneous. 
Such nonaffine deformations reflect the quenched ran- 
domness of the elasticity, and can be derived for a given 
realization of the disorder and a given macroscopic de- 
formation by external stress. Because the deformation is 
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the quantity that is directly measurable in experiments, 
it is useful to derive the relationship between the non- 
affine deformation and the quenched randomness in the 
elastic parameters. Then, by comparing with the RLPM, 
we shall obtain a statistical description of the nonafhne 
deformations, as we shall discuss in Section IV CI 

To study nonafhne deformations, it is convenient to 
take the "state right after linking" [i.e., the state R{z) = 
z] as the reference state, and to re-derive the relaxation 
in the presence of a given deformation A. This is equiv- 
alent to applying the deformation A to the relaxed state, 
and then letting the system further relax for this given 
deformation, as shown in Fig. |5l The relaxed state for 
this given deformation A, which we term the "relaxed de- 
formed state" is described by the deformation z — > z\{z). 
We suppose that 



za{z) = CA • 2 + va{z). 



(4.23) 



For simplicity, we assume that the deformation A is pure 
shear (i.e., det A = 1). 



state right 
after linking 



relaxation 
-t> 



relaxed state 



deformation 



deformed state 
(unrelaxed) 



^ deformation 
I 

^ re-relaxation 



relaxation / relaxed 

' "^-^ I deformed state ; 



FIG. 5: Illustration of the relaxed deformed state. Theo- 
retically, the relcixed deformed state can be reached via two 
routes: in the first, the upper route in this plot, one applies 
the deformation A on the already-relaxed state, and lets the 
system re-relax while keeping the external deformation A to 
the relaxed deformed state; the second route, the lower route 
in this plot, one deforms the system before relaxation is al- 
lowed, and then lets the system relax while maintaining the 
deformation A. These two routes reach the same final state, 
the relaxed deformed state, which characterizes the nonaffine 
deformations that the system undergoes under external de- 
formation. For convenience of calculation, we use the lower 
route to determine the nonafiine deformations. 



and (|E2|) . arriving at 

2{(AaiZi+ {vA)a{z)) J dz2G{z,Z2) 

-2 / dz2G{z,Z2){CAaiZ2,i + {vA)a{z2)) 



-\'oA~^'A~'d,d,{vAh{z)^0. 



(4.24) 



As with the derivation given in Section |IV| we can solve 
this equation pcrturbativcly, to leading order in G*-^' and 
va', see Appendix iGl for details. The result is 



(«a)p 



T 

Pa 



2Dp X'oti\p\^ + 2Dp 



(/a)p, (4.25) 



where and p\, defined in Appendix |Gl are the de- 
formed versions of the projection operators, and ti is 
also defined in Appendix |G] 

In the literature the nonaffine deformations are often 
characterized by the "nonaffine deformation field" wa, 
which is defined in momentum space as 



{wa)p = a ^ • (za)p - Zp , 
= A"^ • {va)p - Vp, 



(4.26) 



where z denotes the relaxed state of the undeformed sys- 
tem (as discussed in Section |IV|, and za is the relaxed 
deformed state. 

Inserting the solution for va into the expression for 
nonaffine deformation field, Eq. (|4.26|) . we have 



iwA)p = 2iC 



{X'a\p\Hi+2Dp)2Dp 

K\p\' ^ 



g 



where 



and 



iK\p\ 



9 



2 + 2Dp)2Dp 



I • • 5p (4.27) 



dp2,a 



P2=0 



P1.P2 ' 



g = ATA 

In the incompressible limit, we have 



(4.28) 



(4.29) 



[wa] 



In Section IV CI we shall compute the mean value and 
disorder correlator of this nonaffine deformation field. 



Next, we use the two stationarity conditions, Eqs. (|4.3j) 
and (|4.7p to solve for the relaxed deformed state. The 
condition (|4.4p for the homogeneous contraction C, is un- 
changed, so we still obtain ^ w 1 — (p/dXo). For the sta- 
bility condition for the random local deformations va we 
follow a similar expansion to the one given in Eqs. (|E1[) 



V. CHARACTERIZING THE PHYSICAL 
ELASTIC QUENCHED DISORDER 

In Section IIVBI we addressed three elastic parameter 
fields — the residual stress cr and the Lame coefficients fip 
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and Ap — which characterize the elastic energy relative to 
the relaxed state and are therefore the physically relevant 
parameters for describing spatially heterogeneous elastic- 
ity. There, we showed how these parameters are deter- 
mined, for a given realization of the quenched disorder, 
by G, i.e., the random nonlocal kernel of the phenomeno- 
logical model, giving the connection in Eqs. (|4.22p . In 
Section Hn] we obtained a statistical characterization of 
G via a comparison with the semi-microscopic RLPM, 
giving the results in Eqs. p.l6p . Thus, we have the in- 
gredients for constructing a statistical characterization of 
the physical position-dependent elastic parameters, as we 
do in the present section. 



A. Disorder averages of the elastic parameters 

Our first step is to determine the disorder average of 
the nonlocal kernel in the relaxed state G. To do this, 
we note that G is related to G via Eq. (|4.17a|) : the 
leading-order expansion of this relationship is given in 
Eq. ()4.19|) . By taking the disorder average on both sides 
of Eq. (I4.19|) . we find that only the first term on the RHS 
survives, because all other terms arc linear in the fiuctua- 
tion part of G and this vanishes upon disorder-averaging. 
Thus, we find that the disorder average of G is given by 

[G(Z1, Z2)] = [G(Z1, Z2)] = Z2), (5.1) 



cragc is given by 

= nokBTd, (5-3) 

with 9 g iven in Eq. ()2.44p . This has been obtained in 
Ref. [15|. This mean shear modulus is linear in tem- 
perature T, reficcting its cntropic nature, a result that 
confirms this aspect of the classical theory of rubber elas- 
ticity. As for the disorder average of bulk modulus A, it 
is obtained via Eq. p.l6c[) . which gives 

IHz)] - ly^nl . (5.4) 

As one might expect, the mean bulk modulus depends on 
the particle number density hq and the strength of the 
excluded- volume interaction i/^. The disorder average 
of the these three elastic parameters of the local form 
of the clastic energy relative to the relaxed state (viz., 
[o'y (z)], [fJ.{z)] and [A(z)]) are all spatially homogeneous 
and isotropic; this is also a consequence of the macro- 
scopic translational and rotational invariance of the ran- 
dom solid state discussed in Section HIE II 



B. Disorder correlators of the elastic parameters 

1. Disorder correlator of the nonlocal kernel 



which means that the disorder average of G is the same 
as the disorder average of G, where we have dropped the 
tilde on z because now we discuss elastic parameters in 
the relaxed reference state only. It is worth noting that, 
as expected, because Ki{zi, Z2) is independent of the cen- 
ter of mass coordinate (zi -I- Z2)/2, the disorder average 
of the nonlocal kernel [G(2;i,Z2)] is translationally and 
rotationally invariant, depending only on |zi — 22]. This 
is a consequence of the macroscopic translational and ro- 
tational invariance of the random solid state discussed in 
Section UTeH 

Second, we determine the disorder averages of the 
position-dependent elastic parameters in the local form 
of the elastic energy relative to the relaxed state, includ- 
ing the residual stress <t, the shear modulus /i, and the 
bulk modulus A. For any given realization of the disor- 
der, these elastic parameters are related to G and Aq via 
Eqs. (|4.22p . Thus, as with the nonlocal kernel, we obtain 
the disorder averages of these elastic parameters via the 
statistics of G. 

The disorder average of the residual stress cr is straight- 
forwardly seen to vanish: 

[a,,{z)]=0. (5.2) 

Thus, the residual stress is a quenched random field with 
zero mean. As for the shear modulus fi, its disorder av- 



Thc nonlocal kernel G characterizes the quenched ran- 
dom nonlocal interactions in the relaxed state. Its statis- 
tics can be described via its moments. In Section IV A[ 
we already determined the disorder average of G; in the 
present section we determine the disorder correlator of 
G. 

To do this, we use Eq. (|4.19p . which relates G to any 
given random configuration of G. Using the disorder cor- 
relator of G, Eq. (|3.16bp . we then arrive at the disorder 
correlator of G, viz., [G{zi, Z2)G(z3, Z4)]. This is a com- 
bination of Gaussian and delta-function factors in the 
separations of the six pairs formed by the four points 
{zi, Z2, Z3, Z4}. The derivation and the momentum-space 
expression of the result for this disorder correlator is 
given in Appendix IH II 

To reveal the universal^ characteristics of the dis- 
order correlator [G(zi, Z2)G(z3, Z4)], we investigate its 
large-distance behavior. The nonlocal kernel itself de- 
scribes a short-distance attractive interaction, because 
the disorder average [G(zi,Z2)] is a short-ranged func- 
tion in |zi — Z2I characterized by the typical localization 
length. Thus, to extract the long-distance behavior of 
[G(zi, Z2)G(z3, Z4)], we take the limit that the two pairs 
{zi, Z2} and {z3, Z4} are far apart from one another, but 
zi is near Z2, and Z3 is near Z4. This is precisely the 
construction of the local description of elasticity of the 
relaxed state introduced in Section IIV Bl featuring the 
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quenched random residual stress and the Lame coeffi- 
cients fields. We shall now discuss the disorder correla- 
tors of these elastic parameters in the local description 
of elasticity. 



2. Disorder correlators of the elastic parameters in the local 
form of the elastic energy 

The elastic parameters in the local form of the elastic 
energy, including the residual stress <t, the shear modu- 
lus fi, and the bulk modulus A, are related to any given 
configuration G via Eq. (|4.22p . Using these relations and 
the disorder corr elator of G, Eq. (fSleb]) . we arrive at 
the disorder correlators of the elastic parameters. The 
details of this calculation are given in Appendix IH 2[ we 
summarize the results in Table U] 
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TABLE I: Long-wavelength variances and covariances of 
the elastic properties of soft random solids in the relaxed 
state. The entry in row R and column C, when multiplied 
Sp+p' Oi yields the connected disorder correla- 
tor [R(p)C(p')]c = [R(p)C(p')] - [C(p')]- 



The correlation function [crcr] features the tensor Aijki , 
which is defined as 



A. 



i'jkl 



T T 



T T 
' PilPjk 



(5.5) 



where the projection operator p''^ is defined in Section lTVl 
The stability condition on the residual stress field cr 
requires that its Fourier transform vanishes when con- 
tracted with the momentum p. It is straightforward to 
see that this feature is obeyed by the correlation function 
[crcr] given in Table HI owing to the structure of A. 

The parameters 6 and 7, on which the correlators in 
Table U depend, are given by 



7 



1 



^fQ - 1 



(5.6a) 
(5.6b) 



The dependence of 9 and 7 on the density of link 
is shown in Fig. [51 The asymptotic behaviors of 9 and 7 
are as follows: 



for 



1(^2 _ 1)3 

for 



for 
for 



> 1; 

if > 1; 



> 1, 



(5.7) 



(5.8) 
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FIG. 6: Plot of Q, Q and 7 as functions of the links density 
parameter r^ . 



where rf is equal to the mean coordination number of 
the particles. Although the connected disorder correla- 
tors of the elastic parameters increase with the density of 
links rf ^ it is worth noting that the relative fluctuations 
are decreasing functions of 77^. For example, the relative 
fluctuation in the shear modulus, defined as [/i/i]c/[/i]^, 
scales as 



i_ 

fl2' 



(5.9) 



which is shown in Fig. [T] This is a decreasing function 
of ?7^, which means that the relative fluctuations of the 
shear modulus actually decrease as links are added. 
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FIG. 7: Plot of the function 0(r?^) = which characterizes 
the relative fluctuations of both the shear and bulk moduli, 
and also the connected correlator of the nonafBne deforma- 
tions, as will be shown in Section TV CI 

It is interesting to look at the real-space behavior of 
disorder correlators of the elastic parameters in the lo- 
cal form of elastic energy. First, it is easy to see that the 
disorder correlators [^(0) ^(r)], [A(0) A(r)] and [^{0) A(r)] 
are short-ranged in real space: more precisely, they are 
proportional to Ss{r), i.e., to a Dirac delta- function that 
has been smoothed on the scale of the short-distance cut- 
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off. This cutoff should be taken to be the typical localiza- 
tion length, in order to validate the Goldstone-fluctuation 
framework for elastic deformations, because the Gold- 
stone fluctuations in the RLPM are long- wavelength, low- 
energy excitations of the random solid state, and these 
do not touch lengthscales shorter than the typical local- 
ization length. 

By contrast, entities involving the residual stress have 
a more interesting spatial correlations: in three dimen- 
sions and at large lengthscales wc find that 



VT^j{Q) akl{r)\c= B^jkl 



(5.10a) 



[a., (0) A.(f)], ^ - ^ ^ (r) ~pI (f)j ,(5.10b) 

[a., (0) A(f)]. = ^^^f^|J^ {p% ir) -Vl (r)) , (5.10c) 

where p^(r) and j3^(r) are, respectively, longitudinal and 
transverse projection operators in real space, which are 
given by 



(5.11) 



The tensor has a complicated structure compris- 

ing terms built from projection operators of r, together 
with various index combinations, and also depends on 
the large-momentum cutoff, which can be identified with 
the inverse of the typical localization length [i^ . 

As the above results show, the mean shear modulus 
and the long-distance behavior of the disorder correlators 
depend only on the link density 77^, the temperature T , 
and the particle density no, and do not depend on the 
details of the link potential. This verifies the argument 
that the shear rigidity of the RLPM is a result of the 
network entropy. 



C. Statistics of nonafRne deformations 

In this section we develop a statistical characteriza- 
tion of the nonafRne deformations of the soft random 
solid state. In Section HV CI we discussed why soft ran- 
dom solids undergo nonaffine deformations in presence 
of a given shear deformation A, and explained how to 
characterize these deformations in terms of the nonaffine 
deformation field wa. 

The nonaffine deformation field is related to any given 
nonlocal random kernel G via Eq. (|4.27p . It is straight- 
forward to sec that the disorder average of the nonaffine 
deformation field vanishes, i.e., [wa] = 0: it is propor- 
tional to the fluctuation part of the quenched random 
nonlocal kernel G^^-* . 

Next, we calculate the disorder correlator of the non- 
affine deformations. For convenience, we take the incom- 
pressible limit, i.e., Aq — )■ 00, in which limit the non- 
affine deformation field wa is given by Eq. (j4.30p . Using 



Eq. (j4.30p . as well as the disorder correlator of the nonlo- 
cal kernel [GG]c, the disorder correlator of the nonaffine 
deformation field wa is found to be 



1 ^ 

no 02 



where 



ti EE Tr(p^g-i), 
t2 = Tr(pLg-ig-i), 
g ^ A^A. 



(5.13a) 
(5.13b) 
(5.13c) 



The dependence of the connected disorder-correlator of 
the nonaffine deformation field [wawa]c on the density 
of links comes through the factor = 7/(0^)^ which is 
shown in Fig.[71 It is evident that, as the density of links 
increases, the system has smaller relative fluctuations of 
its elasticity, i.e., the relative fluctuations in the clas- 
tic moduli decrease, and thus the nonaffine deformations 
also decrease, corresponding to the system becoming less 
heterogeneous. 

The disorder correlator of the nonaffine deformation 
field, Eq. (j5.12p . is consistent with the disorder correla- 
tor of the nonaffine deformations given in Ref. [Tsj . In 
Eq. (3.22) of Ref. [11], the disorder correlator of the non- 
affine deformation u' (which corresponds to ua in our no- 
tation) was found to depend on the random local elastic 
modulus Kiji^i in momentum-space as 



K(9)u'(-9)](x2^^, 



(5.14) 



where 7 represents the appropriate components of the 
tensorial externally applied deformation (i.e., A in our 
notation), represents components of the variance of 
the elastic-modulus tensor, and K represents components 
of the ave rag e of the elastic-modulus tensor. Consistency 
with Ref. [ij] is revealed by taking Eq. (|5.12p . and recall- 
ing that [/x] ~ d and [/i/i]c 7, and therefore that our 
disorder correlator of the nonaffine deformation field can 
be written as 



[(wa)p(wa)-p]c oc 



1 \m\c 



(5.15) 



which exhibits the same dependence on the mean and 
variance of the quen ched random elastic modulus as 
Eq. (3.22) of Ref. [ll does. 

By transforming the disorder correlator of the non- 
affine deformation field (|5.12p back to real space, wc find 
that the large-distance behavior of the disorder correla- 
tor of the nonaffine deformation field [u'a(O) ■ wi\{r)\c is 



proportional to 
long-ranged. 



in three dimensions, which is also 



VI. CONCLUDING REMARKS 

The heterogeneous elasticity of soft random solids 
has been investigated via a semi-microscopic approach. 
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By starting with the Randomly Linked Particle Model 
(RLPM), which describes networks of particles randomly 
connected by soft links, and applying the concepts and 
techniques of vulcanization theory, we have established 
a field-theoretic description of the liquid-to-random-solid 
transition, and have analyzed the corresponding pattern 
of spontaneous symmetry breaking and the structure of 
the associated Goldstone fluctuations. We have identified 
these Goldstone fluctuations as being related to shear de- 
formations of the random solid state and, via this iden- 
tification, we have obtained a statistical characterization 
of the quenched randomness exhibited by the heteroge- 
neous elasticity of soft random solids, which features a 
random nonlocal kernel describing attractive interactions 
between mass-points. 

The heterogeneous elasticity studied via the Goldstone 
fluctuations in the RLPM is a description of the elastic 
properties of the state right after linking (i.e., an elastic 
free energy that takes the state right after linking as its 
elastic reference state). We have shown that, after link- 
ing, the system relaxes to a stable state for any given 
realization of disorder (i.e., for any given heterogeneous 
configuration of the elastic parameters in the state right 
after linking), and this relaxed state, which is a state of 
mechanical equilibrium, is actually the state of experi- 
mental relevance. By solving for the relaxed state for 
any given realization of disorder, and expanding the elas- 
tic free energy for deformations relative to this relaxed 
state, we have obtained an elastic free energy relative to 
the relaxed state (i.e., taking the relaxed state as the new 
elastic reference state). The statistics of the quenched 
randomness in this elastic free energy is subsequently de- 
termined. 

The first statistical moments of the quenched random 
elastic parameters (i.e., the disorder averages of the elas- 
tic parameters), unveil the basic homogeneous macro- 
scopic properties of the heterogeneous elastic medium. 
We have found that the disorder average of the nonlocal 
kernel of attractive interactions is characterized by the 
typical localization lengthscale of the RLPM, which is a 
scale smaller than the lengthscale of the elastic deforma- 
tions that we are considering. Thus, it is reasonable to 
make a local expansion of the elastic energy, relative to 
the relaxed state. The resulting local form of the elas- 
tic energy is a version of Lagrangian elasticity, featuring 
heterogeneous (i.e., spatially randomly varying) residual 
stress and Lame coefficients. The disorder average of 
the residual stress vanishes. The disorder average of the 
shear modulus is found to be proportional to tempera- 
ture, reflecting the entropic nature of the shear rigidity 
of soft random solids. The disorder average of the bulk 
modulus depends on the particle number-density and the 
strength of the excluded- volume interaction. In particu- 
lar, the disorder averages of these elastic parameters of 
the relaxed state are all translationally and rotationally 
invariant, reflecting the macroscopic translational and ro- 
tational invariance of the soft random solid state. 

The second statistical moments of the quenched ran- 



dom elastic parameters (i.e., the spatial correlations of 
these elastic parameters) characterize the fluctuations of 
the quenched randomness in the elastic properties. The 
disorder correlators of the elastic parameters that appear 
in the local form of the elastic energy (relative to the re- 
laxed state) exhibit interesting universal behaviors. In 
particular, the disorder-correlators involving the residual 
stress are found to be long-ranged and governed by a 
universal parameter that also determines the mean shear 
modulus, but the disorder-correlators of the shear and 
bulk moduli are found to be short ranged. 

Because of the heterogeneity present in the elasticity of 
soft random solids, upon the application external stress, 
the system responds by adopting a strain field that is 
nonaffinc (i.e., a strain field that is characterized by an 
inhomogeneous deformation gradient). We have also ob- 
tained a statistical description of these nonaffine defor- 
mations. The disorder average of the nonaffine deforma- 
tions vanishes, and their disorder correlator is also found 
to be long ranged. 

So far, we have studied the first two statistical mo- 
ments of the quenched random elastic parameters of soft 
random solids. The entire probability distribution of the 
quenched random elastic parameters can also be explored 
using the formalism presented here, via the RLPM, and 
one can also progress beyond the local limit of the elas- 
ticity theory. 

This approach to the heterogeneous elasticity of soft 
random solids can also be applied to the setting of liquid 
crystal elastomers, in which the constituent polymers of 
the random network possesses (or are capable of exhibit- 
ing) liquid-crystalline order j27l - l29l . [36j . In liquid crys- 
tal elastomers, the strain field is coupled to the liquid- 
crystalline order, and this produces a rich collection of in- 
teresting phenomena, such as spontaneous sample-shape 
deformation upon changes of temperature, anomalously 
soft modes in the elasticity. The interplay of the hetero- 
geneity of the random network and the liquid-crystalline 
order has interesting consequences: e.g., it can give rise 
to a polydomain structure in the liquid crystalline or- 
der. These interesting topics shall be reserved for future 
studies. 
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Appendix A: Disorder average with the 
Deam-Edwards distribution 



In this Appendix we calculate disorder averages 
weighted by the Deam-Edwards distribution, in partic- 
ular, Z\ and Zi+„ of Section HlBI 

First, we calculate the factor Zi, which is defined as 



The summation over the quenched disorder includes 
two steps: a summation over the number of links A/, and 
a summation over all possible ways of making these M 
links, i.e., of assigning the M links to different collections 
of pairs. Thus Eq. (|A1[) can be written as 
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= ^c{V)(e.p{-£^Y.Ai°\\c,~c,\))j . (A2) 



The mean-field approximation for Zi amounts to taking 
the number density of the unlinked liquid to be N/Vq, 
which is similar to the calculation that yields Eq. (|2.14p , 
and hence we arrive at 



Zi = exp (N In Vq 



2VokBT 2 ^ 
Second, we calculate ^i+„, which is defined as 



(A3) 



= E M! Z^iVo)Z^{Vr. (A4) 



The factor Z^{V)^ can be written in terms of replicas as 



n N 



a=0 j=l 



M 



where Hq is the part of the Hamiltonian Hq (i.e., the 
excluded-volume interaction) for replica a, as defined in 
Eq. (1231) . 

Wc define the Hq average for 1 + n replicas as 
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Y[ n iK - I) ' (^5) Using this notation we arrive at 
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Appendix B: Hubbard-Stratonovich Transformation 



The effective Hamiltonian, Eq. (|2.22p . can be ana- 
lyzed via a Hubbard-Stratonovich (HS) transformation — 
a field-theoretic tool that is often applied to strongly 
coupled models to decouple interactions and develop a 
convenient representation in terms of functional inte- 
grals 0,|3p. 

The version of the HS transformation that we use for 
the RLPM can be illustrated via the following simple ex- 
ample. Consider a statistical-mechanical system having 
the following partition function: 



Z{h) ^ dqe 



Z,{e''^^+'^'>)H,Mym 



where H{q) = Hq{<i) — J<1^ — hq is the total Hamiltonian 
for the variable q, with Ho{q) being the leading-order 
term and Jq^ being considered as a perturbation. (Al- 
though it is just a simple quadratic term, we use it to 
illustrate the method.) The factor Zq = J dqe-""'-'^K 
The term hq denotes the coupling to an external field, 
which generates the statistical moments of q via 



/ \ - ^ 



h=0 



\nZ{h). 



(B2) 



The Jq^ term in the exponent can be decoupled using 
the following version of the HS transformation: 
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{i)'" j du^e-^^-\ (B3) 



with 
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hu) 



4J' 



In this form, the partition function is expressed as an 
integral over the variable w, and the quadratic term in 
the original variable q is now decoupled. If fluctuations 
with large q only appear with very small probabilities, 
as governed by Ho{q), one can expand the ln(e^'^"'?)//(j(g) 
term as a power series in q. Thus, one can obtain an 
effective Hamiltonian Hii-o) via the low-order terms in 
Lo, which has the form of a Landau free energy, and is 
convenient to analyze. 

It is evident that the average of lo, taken with the sta- 
tistical weight defined by 7^(0;), equals the average of q, 
taken with the statistical weight defined by H{q): 
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dh 
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\nZ{h) = {q)H{q)- 



(B5) 



Thus, the statistical mechanics of q can be examined by 
studying the statistical mechanics of w. In the Hamilto- 
nian T-L(lu), q appears linearly; therefore, in cases in which 
g is a variable that involves a summation over many par- 
ticles, this method will allow us to decouple the problem 
into a single-particle one, as will be seen in the follow- 
ing application of the HS transformation to vulcanization 
theory. 

In the RLPM, the partition function we are going to 
decouple is Eq. ^JQ: 
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(B4) with 
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The field Qp = (1/A'') ^ '^ '^^ is a complex field, so we apply the following equalities for the complex variables q 

and w: 

e'-'l^l'^^^ Jd{Rcuj)d{lmuj)e~J\'^\'+^'-'^'""^\ (B8a) 
= ^ y"d(Rew)d(Imc^)e-^l'"l'+2^^°9"*, (BSb) 



e+ 



noticing that the product Re((7aj*) = (Re (7) (Re w) + {lniq){lmLd). We use Eq. (jB8a[) for the HS transformation for 
the LRS fields, and Eq. (IB8bp for the HS transformation for the HRS field, and hence arrive at the form (For more 
details, see Section V in [3|.) 

= VQpH Vn^e (B9) 

n 



with 



Hnp;,np] = ^^^^^^ ^i>^-P^p^''^ +^^r!p,of7_p,o +^^X^r!p,.f^_p,»-iVfcBTlnSo, (BIO) 

2F"Aq p^hrs ° P P a=l 



where the iVfesTlnSo term is analog of the ln(e^'^"'')i/^(g) term in Eq. (|B4p and, using Q_p = (1/^) e'-^'^^ and 

Q^p - (1/7V) Ef=i e'^-^^° , we have 
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I 

In this form it is evident that the N particles are actually point order parameter (|2.40p into the Hamiltonian (|2.3ip . 
decoupled. Notice that in Eq. (jB9p the functional inte- 



grals Jvrip na=o have carefully chosen prefactors 
[as in Eq. (jB8|) ] , to ensure that the integration is properly 
normalized. 



Appendix C: Hamiltonian of the stationary point 



In this Appendix we calculate the value of the Hamil- 
tonian at the stationary point by inserting the stationary 



The first term in the Hamiltonian is 
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_ jV7?^fcgrQ^(A(P))" jV7?^fcBTg^(AW)" /■ f / 1 1 

2T/« 2F" ^ y^.r^^ ^2ri 2t2 2kBT)i ' ^ 

where = {z, ^z, . . . , Cz}. The sum 'Ylip^HRS ^^"^ changed into the because the order parameter we have 
inserted vanishes for p gLRS. We have also changed momentum summation into an integral by using y^yn X]j3 = 

J (27r)(i + ")<i • ______ 

The free energy of the system is related to the 0(n) term of this Hamiltonian, as given by Eqs. (|2.10[ [^?TT|) . Thus, 
we make the small-n expansion. It is straightforward to sec that the 0(1) terms cancel, and that the leading-order 
term is given by 0{n) 

^0 pGHRS 

Nri^keTQ'^ ( d,, ,^ , d [ /I 1 



Similarly, we can calculate the In^o term. By inserting the saddle-point value of the order parameter into Sg, and 
summing (or, more precisely, integrating) over momentum p, we have 



= e-'''0(^c'"/vry'rfgexp{r;2g(A(°))" J ^ (|-) "^g-J }, (C3) 



where f = {^^ + a^) ^ . We then Taylor-expand the exponential (keeping all orders) and integrate out c to get 
So = e'^'^i^^!"fyrv,V-{l + ,^Q(A(°)/F)" + E ^tl / - • - '^'» 
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n/T,-\^T— / 27r "i"2(^c,i-'c,2) +■■■ 1 

, ,„,^,(2i) (jmw) ' }■ <'^'" 

where in the exponential the terms following (z^^i — z^,2)^ include all pairs of the m variables [there are m(m — l)/2 
such terms]. By using {z(^,i — Z(^^2)^ = {1 + n('^){zi — z-if' [recall that z^^i is a (1 -|- n)d-dimensional vector, and that 
zi is a d-dimcnsional vector] , the integration J dz\ ■ ■ ■ dzm can be readily performed, and we thus obtain 
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From this, — In can obtained by making the small-n expansion, using the following equality: 
ln(x + ny + 0{n^j) = ln(x-(l + n{y/x) + 0{n'^))) = Inx + n{y/x) + 0{n'^), 

so that we have 

-In So = -lnyo+n{-lnl/+(77'(3 + e-''''5-l)(^(ln(2^) + C')-lnl^ 
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Therefore, the small-n expansion of the stationary-point Hamiltonian H^^"^ to 0{n) is given by 
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where the parameter 6 is defined via 
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Appendix D: Hamiltonian of the Goldstone 
deformed order parameter 

In this Appendix we calculate the value of the Hamil- 
tonian for the Goldstone-deformed order parameter by 



inserting the Goldstone-deformed order parameter (j2.59p 
into the Hamiltonian (|2.3ip . following a calculation sim- 
ilar to that in Appendix [C] To obtain a description 
of the elasticity, we shall expand the Hamiltonian for 
small deformations, specifically in a series in the small, 
scalar variable that characterizes the replicated deforma- 
tion field, viz., *fl(zi,Z2) = {R{zi) - R{z2)f - (1 + 



The quadratic term in the Hamiltonian is given by 
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Next, we expand for small 'i'R, adopting the notation Z2) = {R{zi) — R{z2)Y ^ (1 + — ^2)^- Note that 

^ u is not related to the deformation relative to the stationary point, this stationary point being characterized by 
the mean positions of the replicas of the particle z^ = {z, C,z, . . . , C,z}. Instead, ^ describes deformations relative to 
the "state right after linking" (i.e., prior to relaxation), this state being characterized by the mean positions of the 
replicas of the particle z = {z, z, . . . , z}, as discussed in Sections III E 31 and IIVI The expansion of the quadratic term 
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for small is given by 
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The small-n expansion on this quadratic term is then given by 
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The calculation for the In So term is similar to the above calculation of the quadratic term. The expansion in small 
quantity r reads 
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where the summations that wc have abbreviated with • • • include all pairs formed by {zi, . . . , z,„}. Next, we expand 
for small n, keep terms to 0{ri) in the In So term, assuming that both 4*/^ and contain 0{n) terms. After a 
tedious calculation we have 
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To further simplify the expression, first consider the r) terms in the expansion, ~ '^^'2^i + ^-^-+^'\" ■ ^^'^ ^''^t term 
has a factor of ii{zi, Z2) = (-R(^i) — -^(^2))^ ~ (1 + ''t-){zi ~ ^2)'^, which only involves two variables zi and Z2, so we 
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total, there are such terms (i.e., the number of pairs among m variables). Thus, the 0{'^ii) term in — InSg 



can integrate out the other (m — 2) variables, i.e., 23, ... , Zm- (Of course, for m = 2, no integrals are needed.) In 

total, there 
is given by 
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where in the last line here we have used the fact that {zi, Z2, ■ ■ ■ , Zm} appear symmetrically, so that the "^^"^~^'> terms 
are identical. 

Similarly, for the 0(^'|j) terms in the expansion in Eq. (jPSp there are terms involving two points, such as ^ij(zi, 2:2)^, 
three points, such as "i!ii{zi, Z2)^r{zi, 23), and four points, such as ^'/{(zi, 22)^^(2^3, (Of course, for m = 3 there 
are no four-point terms, and for m = 2 there are no three or four points terms.) Thus, the 0{^^) terms can be 
written as 



■Tif2 '5^(^1,22) + 
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+m[m - l)(m - 2)fif2 fa^-ifX^i, •Z2)*fl(z2, 23) 
m{m — l)(m — 2)(m — 3) _ _ - _ 
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Following a calculation similar to that in Eq. (|D6|) . we can integrate out the integration variables that are not present 
in V^ij, and thus obtain the 0(4'^) term in IuSq. 

Summing up the contributions from the quadratic term and the InSg term, we arrive at the Hamiltonian of the 
Goldstonc deformed state: 

H^f^^H^i^^+Hl-, (D8) 

with being the Hamiltonian of the stationary point, and the increase of the Hamiltonian due to Goldstone 

deformation is given by 

if*" = -nNkeT^iC^ -1) + ^ j dzidz2Ki{zi,Z2)^B,{zi,Z2) 

- „ / ^ [ dzidz2dZ3dZ4K2{zi,Z2,Z3,Z4)'iii{zi,Z2)'iR{z3,Z4). (D9) 

okb-l j 

The first term here, i.e., —NkBT^{('^ — 1), is present due to the fact that the expansion variable \E'ij;(zi, Z2) measures 
departures from the state right after linking, not the stationary point, as we have previously discussed. Note that -ff*" 
involves only the energy of shear deformation, because the Goldstone modes contain only pure shear deformation. 
The energy of volume variations is in the stationary-point Hamiltonian part, which contains the variable contraction 
parameter C. 
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The kernels in Eq. (|D9|) are given by 
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Strictly speaking, the kernel A'2 should be symmetric under the exchanges of variables zi 4-> Z2 or Z3 4-!> Z4. Here, 
to save space, we have written the above non-symmetric form. The true (i.e., symmetric) form can be recovered by 
averaging: 

i^2 (21,2:2,23,24) ^(^2(21, 22, 23,24) +A'2(2i, 22,24,23) + Ar2(22,2i,23,24) +1^2(22,21,24,23)). (D12) 



Appendix E: Relaxation of the phenomenological 
elastic free energy for a given realization of disorder 



Using this expansion, we have 



In this Appendix we solve the stationarity condition 
for the random local deformations v. First, we need to 
calculate the variation of the "bulk term" which can be 
expanded, to leading order in small u, as (46j 
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(^) ^ -( 



det(^C% + djv^{z)^ 
C''(l + r'%'.(2)). (El) 



+e-'-^d,v,i^)djVj{z). 



Thus, the stationarity equation reads 



(E2) 
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= 2{(:Za + Va{z)) / dZ2G{z,Z2)-2 dZ2G{z , Z2){C,Z2,a + Va{z2)) - Kda{d.m{z)) , 



where 



-^0 = ^oC 
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We take the disorder average of the nonlocal kernel G*^°^ to be a zeroth-order quantity, and the fluctuation part G*-^^ 
to be a first-order quantity and. thus, v{z) is also of first order. The he zeroth-order equation is then 



Q^2Za dZ2G^°\z -Z2)-2 / dZ2G^°\z - Z2)z2,a 
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which is automatically satisfied, given that G^'^\z — Z2) is even in (z — Z2). 
The first order equation reads 
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We address this equation in momentum space. Wc define the following Fourier transforms (on a specific finite 
volume — the volume of the state right after linking, viz., Vq): 



G^°'> = j dxe-'P^G^°\x), 
so that the momentum-space stationarity equation becomes 
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Strictly speaking, the derivatives here should instead be 
understood as difference quotients because we are using 
a finite-volume version of the Fourier transform; but for 
convenience we write it a derivatives. 

Equation (|E8p can be written in the tensorial form 



{2(^)1 + A;,pL} . \p\^v{p) ^ /(p), (E9) 



where 
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This quantity fa,pi is actually the random force in the 
state that is contracted but has not yet been equilibrated 
for the randomness, and 
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Furthermore, I is the d-dimensional identity matrix, and 
the projection operators in momentum space, pfj and pjj, 
are defined as 



Pii = {PiP])/P^ 



(E12) 



They satisfy the following relations: 
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(pT)2^p^ p^.p^^O. (E13) 



In the following we shall use bold-face letters to denote 
rank-2 tensors, and letters with an overhead arrow (such 
as v{p)) to denote a vector, when needed. 

By this decomposition we arrive at the solution to 
Eq. (lE9l) : 



P'^ ■ fp P^ ■ fp 
2Dp X'o + 2Dp " 



(E14) 



Notice that the second term is much smaller than the 
first term, due to the large bulk modulus Aq. In the 
incompressible limit (i.e., Aq 00), we have 



fp 



2Dr, 



(E15) 



which is a purely transverse field, meaning that it satisfies 
Pi ^i,p = or, equivalently, di Viix) = 0, which is the only 
deformation allowed in an incompressible medium. 
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Appendix F: Re-expanding the elastic energy 
around the equilibrium reference state 



reference space labeled by x, and as a scalar in the target 
space labeled by R. 



In this Appendix we re-expand the elastic energy for 
deformations relative to the relaxed state, z = C,z -\- v{z)^ 
as discussed in Section flVBI The small variable in this 
expansion is the deformation field u{z) relative to the 
relaxed state. Furthermore, to obtain a continuum de- 
scription of the elasticity, we adopt a notation involving 
the strain tensor eij[x): 

= l{d.u,ix)+dM-)+dM-)dM-))- (Fl) 

where ^ij{x) = dRi(x) / dxj is the deformation gradient 
tensor. This strain tensor transforms as a tensor in the 
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1. Expanding the nonlocal kernel G in the relaxed 
state 



The definition of G, given in Section IIVB| is 



G(zi,Z2) = G{z{zi),z{z2)). 



(F2) 



It can be expanded for small v to yield a direct expression 
for G. In momentum-space this reads 
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dzidz2J{zi)J{z2)e'''P''-<''+''^''^')-'P'^<''+''^''^^G{zi,Z2), 



(F3) 



where, in the first line, z(zi) is the mapping of a mass point zi in the relaxed state back to the position z(zi), at 
which it was located in the state right after linking. Inserting in the expressions for C and v{z), given in Eqs. (|4.5[ 
14. 9p . and keeping terms to O((l/Ao)°) [which gives ^ ~ 1 and J^{z) ~ 1], we can expand v down from the exponent, 
and keep terms to first order in G*-^' (noting that v is the same order as G^^-*), and thus arrive at 



dzidz2{l ~ ipiv{zi)-ip2v{z2))e-'P^''-'P-^^{G^''\zi,Z2) + G^^\zi,Z2)) 



Gplp. + GpL - ' I dz,dz2{pMzi)+P2v{z2))e-^P^^^-^i'^^^G^°\z,,Z2) 



'Gpi,P2 + 'Gi>UP2 - • '"ipi+P2)Gp^^ +P2 ■ Vipi+P2)Gp, 



(o)^ 



(F4) 



2. Local expansion of the harmonic attraction 



In this section we make a local expansion of the nonlocal term in the elastic free energy of the equilibrium reference 
state, Eq. (|4.16p . i.e., the term 



Tnonlocal =^ j dzidz2J{zi) ^ J {Z2) ^ G{zi, Z2){\^R{zi) R-ih)]^ - \zi ^ 

For convenience of notation, we define the following change of variables: 

z = Zi, 

V = Z2- Zi, 

M{z,y) = G(zi,Z2), 

so that the nonlocal kernel in the relaxed state, Eq. (|4.19p . can be written (in momentum space) as 

with the definition of gp.q, and then its leading-order expansion in momentum, being given by 

i{Gf-Gf_^) tip^-2p-q) 



9p,< 



2(4°^ -Gf) 



2p 



n2 



(F5) 



(F6) 
(F7) 

(F8) 
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The local expansion of Eq. (|F5[) then becomes 

Tnoniocai = ^ J dMij M {z , y) (jR{z) ~ R{z + y)\^ - \y\^^ 



~ -y dS{d,RiiS)djRi{S)-S^j) J dyy,yjM{S,y), (F9) 

where the factor J^{zi)~-^J'{z2)''^ is ignored because its difference from unity is of O{l/Xo). Now it is straightforward 
to express the elastic energy Fnoniocai in the standard form of Lagrangian elasticity using the strain tensor eij{z) = 
^{dMi5)djRiiz)-S,j). 

The complete expression of the local form of the elastic energy for deformations relative to the relaxed state, 
including the contribution from the bulk term, will be calculated in Appendix IF 41 



3. Expansion of the bulk term 

The "bulk term" in the elastic free energy Eq. ()4.16|) is given by 

•^0 f ^,..jdR^{S) 



Tbulk - 2 



'-J 

The determinant in this equation can be expanded using the strain tensor e: 



d~zj{z)-^[j{z) dct(^||;^) - l} . (FIO) 



det(^l^) = det(A(z)) = {det(I + 2e(i))}^/' = e^Trin(i+2.(i)) 



1 + Tre(z) - Tri{zf + i(Tre(5))^ + 0{i{S)^). (Fll) 



Thus, we have 



^0 /"j-^/ T/ j9Ri{z) 

t bulk — ~ - - " 



dSJ{z)-'{j{z)dct[^^) -l} 

~ d2J{z)-'[{Jiz) - If + 2{J{z) - l)J(z)Tre(z) 

- 2{J{z) - l)J(z)Tre(f)2 + (2J(z) - 1) J(z)(Tre(5))2|. (F12) 
Inserting the solutions for C and v, given in Eqs. (|4.5[ 14.9^ . into the Jacobian ^{z) = we arrive at 



Tbuik - / dz{Tr(<T'(z) • e(z)) + Ai(z)Tre(z)2 + ^(Tre(2))2}, (f13) 



with the elastic parameters (in momentum space) being given by 

■ip- f, 



t^p = pVoSp - (F14b) 
Ap = XoVo6p + 2{'-^^-pVo5py (F14c) 

4. Local form of the elastic energy relative to the relaxed state 

Summing up the contributions from the nonlocal term F„oniocai and the bulk term Fbuik to the elastic free en- 
ergy (j4.16p . we arrive at the local form of the elastic energy for deformations relative to the relaxed state: 

F = 1 dz{Tr(<T(f) • i{S)) + /z(z)Tre(5)2 + ^(Tre(5))2}, (F15) 
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with the clastic parameters being given by 

Gp-q,p + —r^ ~ —r^ {P^PJa.p + PoP^a,p) ' (Fl^a) 



fip^pVoSp^'-^, (F16b) 

Ap = Aoyo-^p + 2{^^ - P^o'Jp}- (F16c) 

I 



Appendix G: Relaxation of the deformed state tion IIV C[ and thus obtain information about nonaffine 

deformations. The stationarity condition is given by 

In this section we solve the stationarity equation with 
a given macroscopic deformation A, as discussed in Sec- 



2{CAa^Z^ + («A)a(z)) J dz2G{z,Z2) ~ 2 j dz2G(z, Z2 ) (CAa»^2,^ + KjaC^a)) - X'^t^^ K~^' drdj{v K)b{z) = 0. (Gl) 

We take G'-^' to be of zcroth order, and G^^^ and wa(z) to be first-order quantities. Thus, the zeroth order equation 
reads 

- 2Ka^Z^ [ dzaG^^^ (z - Z2) - 2A„ / dz2G(°) (Z - Z2)z2,i, (G2) 



which is already satisfied given G^'^^z — Z2) is even in (z — 22)- 
The first-order equation reads 

CAa,Z, f dz2G^^'>{z,Z2) + {vA)a{z) f dz2G^''\z , Z2) - f dZ2G^^\z , Z2)CAa^Z2,^ - I dZ2G^'''> {z , Z2){va) a{z2) 



-^A-a'Aj,'d,d,{vA),{z) = 0. (G3) 
We can solve this equation in momentum space, where it is expressed as 



= jCA^-^G^o - ^CAa^J^ Gl\l^ + (G[,°^ - G^^^^)ivA)a,p, + ^A,-iA7,Vi,,pi,,«6,p,. (G4) 
Opi i (-'P2A P2— U ^ 



Writing this equation in tensorial form, we have 



2(^)1 + a;,a-Va-i ■ bp(«A)p = (/a)p, (G5) 



where g = A-^'A is the metric tensor, and 

Dp ^ g[,°)-GW, (G6a) 
iU)a,p, ^ -2c[^Aa^-^Gl^l,~^Aa^^Gl%,). (G6b) 

I 



To solve Eq. (|G5|) . it is useful to define the A-deformed It is straightforward to verify that they obey 
versions of the projection operators, i.e., 

1 (^7.) (Pa)'=PA' (pD'^pL P^-pI-O. (G8) 

" Tr(pLg-i) P ' > 

pI ^ I-Pk- (G7b) 
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By using these projection operators we can write 
Eq. ((U5|) as 

{^pI + + Kh)p\}- \p\'Mp = (/a)p,(G9) 

where we have defined 

h = Tr(pLg-i). (GIO) 
Thus, it is straightforward to arrive at the solution: 

I 



Appendix H: Correlation functions of the elastic 
parameters in the equilibrium reference state 



1. Correlation function of the non-local kernel in 
the equilibrium reference state 



The nonlocal kernel in the equilibrium reference state 
G is related to the nonlocal kernel in the state right after 
linking G via Eq. (|4.17aP : to leading-order in the small 
quantity v we have 



.(1) 



(0) 



(HI) 



By using this relation, we can derive the correlation function of G from the correlation function of G [which is given 
m Eq. ()3.16bp ]. and thus we arrive at the correlation function 



[Mp,,q,Mp,^g,], 

2kBT 



1 1 

n T2 



1 ( 



/2 



N g (yy^Q)- m(m - 1) 



N 



m— 2 



Ti , . ...Tr, 



TlT2 



2 / (ti +T2 +92 I 

g 2x1 ■f2 



(Tl+x2)|pi--gi+g2r 
. g 2x11-2 



^ Yin] 



m— 3 

1 ( l<;il^ I lll+il2l^ I l92 



T2 ' X3 / g 2 V XI X2 X3 



,...,T^ (n H h 



_1 ( |pi-qil I |pi-gi + 92l I Iggl 

g 2 \ XI X2 T3 



I |pi-gi+g2l I |pi+g2l 



>)/4 



-e — V 

(fcBT)2 m! 4 



-,^Q (7?^Q)'"m(m-l)(m-2)(m-3) 



m=4 



2 V Tl ?2 ?4 



T r ^rQ /■ / 1 1 

L 2kBT J \Ti T2 



+ 



+ 



X 




^-pim 








X 


^pi,qi 






{kBTf 


X 


^pi.gi ( 



V2 



^ ml 9 



r„ (Tl + ---+r„)2 

(xi+T2)|9l| 



(T-l-|-x2)|g2l (T-i+x2)|pi+g2r , , 

> ; — m(m— l)(m — 2) 



m— 3 



(ti +-r2)l -Pl+gi I" , 

e 2Tif2 )/2 



(n H h T„ 



(H2) 
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where we have used the notation M(z,y) = G{zi,Z2) form of the elastic energy for deformations relative to the 
defined in Appendix IF 21 relaxed state. 



2. Disorder correlators of the elastic parameters in 
the local form 

In this appendix we calculate the disorder correlators 
of the quenched random elastic parameters in the local 



First, we calculate the disorder correlator of the resid- 
ual stress [<J<j]c- The residual stress cr in the relaxed state 
is related to the nonlocal kernel G via Eq. (|4.22ap . Thus, 
the correlator of the residual stress can be expressed as 



d 



11,1 
2 

\Pi\ 
2 



d 



11=0 dq: 



92=0 



:{Pl,rVlj{Pl) +PijpL(pi) +Pl,aPybl)) 



d 



dq- 



2.k 



<Jl=0 



92=0 



-■^j,Plfi^l,P2,q2\c 



■^|^^py^(pi,iPaj(Pi) + PijP«(Pi) + Pi,aP5(Pi))(Pi.fcP (H3) 



r 



where we have defined Nj^p^q = dMp^q/dgj, and the no- 
tation AI{z,y) = G(zi,Z2) is defined in Appendix IF 21 

We then insert in the disorder correlator 
[Mp-^.q-^Mp2^q2]c, given in Eq. (|3.16bp (in the form 
of [GG]) into Eq. (jH3p . After a lengthy calculation, and 
making use of the identity 



+771(771 — 1) 



n T2 



1, (H4) 



WG arrive at the correlator 



Spi+P2jj—fy^i'^PijPll + PilPjk + PikPjl), (H5) 

= -\7fQ'^+7fQ+e-'^^^ -I is given in Eq. ((C9|) . 
By following a similar scheme, we have also calculated 
the other disorder correlators and cross-correlators of the 



where 6 



quenched random elastic parameters in the local form 
of elasticity of the relaxed state. Hence, we arrive at 
the correlators of the shear modulus and bulk modulus, 
which read 



[MpiMp2]c =l5p^+p2N{kBTf, 



(H6a) 
(H6b) 



in which the dimensionless scale factor 7 is given by 
3 



^-"^ -l + TfQ + iTi'Qf. (H7) 
We also arrive at the cross-correlators, which are given 

by 



-2N{kBTfe5p,+p,pl{p^), (H8a) 
m{kBTfe5p,+p,^l{pi), (H8b) 
-2N{kBTfl5p,+p,- (H8c) 



the scale factor 6 is defined in Eq. (|C9[I . 
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